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ABSTRACT 



We study Milky Way kinematics using a sample of 18.8 million main-sequence 
stars with r < 20 and proper-motion measurements derived from SDSS and POSS 
astrometry, including ~170,000 stars with radial-velocity measurements from the 
SDSS spectroscopic survey. Distances to stars are determined using a photomet- 
ric parallax relation, covering a distance range from ~100 pc to 10 kpc over a 
quarter of the sky at high Galactic latitudes {\b\ > 20°). We find that in the 
region defined by 1 kpc < Z < 5 kpc and 3 kpc < i? < 13 kpc, the rotational 
velocity for disk stars smoothly decreases, and all three components of the veloc- 
ity dispersion increase, with distance from the Galactic plane. In contrast, the 
velocity ellipsoid for halo stars is aligned with a spherical coordinate system and 
appears to be spatially invariant within the probed volume. The velocity distri- 
bution of nearby (Z < 1 kpc) K/M stars is complex, and cannot be described 
by a standard Schwarzschild ellipsoid. For stars in a distance-limited subsam- 
ple of stars (<100 pc), we detect a multimodal velocity distribution consistent 
with that seen by HIPPARCOS. This strong non-Gaussianity significantly af- 
fects the measurements of the velocity ellipsoid tilt and vertex deviation when 
using the Schwarzschild approximation. We develop and test a simple descriptive 
model for the overall kinematic behavior that captures these features over most 
of the probed volume, and can be used to search for substructure in kinematic 
and metallicity space. We use this model to predict further improvements in 
kinematic mapping of the Galaxy expected from Gaia and LSST. 

Subject headings: methods: data analysis — stars: statistics — Galaxy: disk, halo, 
kinematics and dynamics, stellar content, structure 
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INTRODUCTION 



The Milky Way is a complex and dynamic structure that is constantly being shaped 
by the infall of matter from the Local Group and mergers with neighboring galaxies. From 
our vantage point inside the disk of the Milky Way, we have a unique opportunity to 
study an ~ L* spiral galaxy in great detail. By measuring and analyzing the properties 
of large numbers of individual stars, we can map the Milky Way in a nine-dimensional 
space spanned by the three spatial coordinates, three velocity components, and three stellar 
parameters - luminosity, effective temperature, and metallicity. 

In this paper, t he third in a serie s of related studies, we use data obtained by the Sloan 



Digital Sky Survey (lYork et al 



20001 ) to st udy in detai l the d istribution of tens of millions 



of stars in this multi-dimensional space. In 



Juric et a. 



20081. her e after 



Ivezic et al. 



08), we examined 



(l2008al . hereafter 108) we 



the spatial distribution of stars in the Galaxy, and in 
extended our analysis to include the metallicity distribution. In this paper, working with a 
kinematic data set unprecedented in size, we investigate the distribution of stellar velocities. 
Our data include measurements fro m the SDSS astrometri c, photometric, and spectroscopic 



surveys: the SDSS Data Release 7 (lAbazajian et al. 



20091 ) radial-velocity sample includes 



~ 170,000 main-sequence stars, while the proper-motion sample includes 18.8 million stars, 
with about 6.8 million F/G stars for which photometric metallicity estimates are also 
available. These stars sample a distance range from ~ 100 pc to ~ 10 kpc, probing much 



farther from Earth than the 



[e.g., 



Dehnen &: Binney 



3IPPARC0S sample, w hich covers only the nearest ~ 100 pc 



1998 



Nordstrom et al. 



20041 ). With the SDSS data set, we are 



offered for the first time an opportunity to examine in situ the thin/thick disk and disk/halo 
boundaries over a large solid angle, using millions of stars. 

In all three of the papers in this series, we have employed a set of photometric parallax 
relations, enabled by accurate SDSS multi-color measurements, to estimate the distances to 
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main-sequence stars. With these distances, accurate to ~ 10 — 15%, the stellar distribution 
in the multi-dimensional phase space can be mapped and analyzed without any additional 
assumptions. The primary aim of this paper is thus to develop quantitative understanding 
of the large-scale kinematic behavior of the disk and halo stars. From the point of view 
of an observer, the goal is to measure and describe the radial-velocity and proper-motion 
distributions as functions of the position in, for example, the r vs. g — r color-magnitude 
diagram, and as functions of the position of the analyzed sample on the sky. From the 
point of view of a theorist, we seek to directly quantify the behavior of the probability 
distribution function, p{vr,V(I„vz\R, 4>, Z,T, L,[Fe/B.]), where (f_R,f<^,fz) are the three 
velocity components in a cylindrical coordinate system, {R, 0, Z) describe the position of a 
star in the Galaxy, and T, L, and [Fe/H] are its temperature, luminosity, and metallicity, 
respectively ("|" means "given"). 



Th is a different approa c 



model (iRobin &: Creze 



1986 



1 than that taken by the widely-used "Besangon" Galaxy 



Robin et al. 



2003 



and references therein), which attempts to 



generate model stellar distributions from "first principles" (such as an adopted initial mass 
function) and requires dynamical self-consistency. Instead, we simply seek to describe the 
directly observed distributions of kinematic and chemical quantities without imposing any 
additional constraints. If these distributions can be described in terms of simple functions, 
then one can try to understand and model these simple abstractions, rather than the full 
voluminous data set. 

As discussed in detail by JOS and 108, the disk and halo components have spatial 
and metallicity distributions that are well-fit by simple analytic models within the volume 
probed by SDSS (and outside regions with strong substructure, such as the Sgr dwarf 
tidal stream and the Monoceros stream). In this paper, we develop analogous models that 
describe the velocity distributions of disk and halo stars. 
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Schwarzschild 



ipsoid al approximation (a three- 



19791 ) for describing the velocity 



Questions we ask include: 

» What are the limitations of the Schwarzschild el 
dimensional Gaussian distribution 
distributions? 

» Given the increased distance range compared to older data sets, can we detect spatial 
variation of the best-fit Schwarzschild ellipsoid parameters, including its orientation? 

» Does the halo rotate on average? 

» Is the kinematic difference between disk and halo stars as remarkable as the difference 
in their metallicity distributions? 

» Do large spatial substructures, which are also traced in metallicity space, have 
distinctive kinematic behavior? 



Of course, answers to a number of t 



excellent reviews, see 



Gilmore et al. 



l ese questions are kno^vn to some extent (for 



1989 



Majewski 



1993 



Helmi 



2008. for context and 



references, see also the first tw o papers in t 



least since the seminal paper of 



Eggen et al 



th is seri es). For example, it has been known at 
( 119621 ) that high-metallicity disk stars move on 
nearly circular orbits, while many low-metallicity halo stars move on eccentric, randomly 
oriented orbits. However, given the order of magnitude increase in the number of stars 
compared to previous work, larger distance limits, and accurate and diverse measurements 
obtained with the same facility, the previous results (see 108 for a summary of kinematic 
results) can be significantly improved and expanded. 

The main sections of this paper include a description of the data and methodology 
(§2), followed by analysis of the various stellar subsamples. In §3, we begin by analyzing the 
proper-motion sample and determining the dependence of the azimuthal (rotational) and 
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radial-velocity distributions on position for halo and disk subsamples selected along / = 0° 
and / = 180°. The spectroscopic sample is used in §4 to obtain constraints on the behavior 
of the vertical- velocity component, and to measure the velocity-ellipsoid tilts. The resulting 
model is then compared to the full proper-motion sample and radial-velocity samples in §5. 
Finally, in §6, we summarize and discuss our results, including a comparison with prior 
results and other work based on SDSS data. 



2. DATA AND METHODOLOGY 



T 

work 


le characteristics of t 


le SDSS imaging and spectroscopic data relevant to this 


Fukueita et al. 


1996: 


Gunn et al. 


1998: Hoee et al. 


2001: Smit 


1 et al. 


2002: 


Stoue 


iton et al. 2002: 


Pier et al. 2003: I 


vezic et al. 2004: ' 


Fucker et al. 


2006: 


Gunn et al. 


2006: 


Abazaiian et al. 


2009: 


Yannv et al 


20091) are described in detail in the first two 



papers in the series (J08, 108). Here, we only briefly summarize the photometric-parallax 
and photometric-metallicity methods, and then describe the proper-motion data and their 
error analysis. The subsample definitions are described at the end of this section. 



2.1. The Photometric Parallax Method 

The majority of stars in the SDSS imaging catalogs (~ 90%) are on the main sequence 
(J08 and references therein) and, using the broadband colors measured by SDSS, it is 
possible to estimate their absolute magnitude. Briefly, the r-band absolute magnitude, 
Mr, of a star can be estimated from its position on the stellar locus of the Mr vs. 
g — i color-magnitude diagram. The position of this stellar locus is in turn sensitive to 
metallicity, so we must apply an additional correction to the absolute magnitude. A 
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maximum-likelihood implementation of this method was introduced and discussed in detail 
in JOS. The method was further refined by 108, who calibrated its dependence on metallicity 
using globular clusters. 

We estimate absolute magnitudes using equation A7 in 108, which corrects for age 
effects, and equation A2 in the same paper to account for the impact of metallicity. The 
resulting distance range covered by the photometric parallax relation depends upon color 
and metallicity, but spans ~ 100 pc to ~ 10 kpc. Based on an analysis of stars in globular 
clusters, 108 estimate that the probable systematic errors in absolute magnitudes determined 
using these relations are about 0.1 mag, corresponding to 5% s ystematic di s tance errors (in 
addition to the 10 — 15% random distance errors). In addition, ISesar et al.l (120081 . hereafter 
S1J08) used a large sample of candidate wide-binary stars to show that the expected error 
distribution is mildly non-Gaussian, with a root-mean-square (rms) scatter in absolute 
magnitude of ~ 0.3 mag. They also quantified biases in the derived absolute magnitudes 
due to unresolved binary stars. 



2.2. The Photometric Metallicity Method 

Stellar metallicity can significantly affect the position of a star in the color-magnitude 
diagram (there is a shift of ~ 1 mag between the median halo metallicity of [Fe/H] ~ — 1.5 
and the median disk metallicity of [Fe/H] ~ —0.2). SDSS spectroscopy is only available for 
a small fraction of the stars in our sample, so we adopt a photometric metallicity method 
based on SDSS u — g and g — r colors. This relation was originally calibrated by 108 using 
SDSS spectroscopic metallicities. However, the calibration of SD SS spectroscopic metallicity 



chang ed at the high-met allicity end after SDSS Data Release 6 (lAdelman-McCarthy et al. 



20081 ) . Therefore, we recalibrate their expressions as described in the Appendix. The new 



calibration, given in equation (Al), is applicable to F/G stars with 0.2 < g — r < 0.6 
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and has photometric-metallicity errors that approximately follow a Gaussian distribution 



SPSS spectroscopic metallicity (Beers et al 



2008a 



AUende Prieto et al. 



2006 



- dex systematic uncertain 


;ies in 


AUende Prieto et al. 


2006 




Lee et al. 



20081 ) are inherited by the photometric metallicity estimator. 
We emphasize that photometric metallicity estimates are only robust in the range 
—2 < [Fe/H] < (see Appendix for details). 

For stars with g — r > 0.6, we assume a constant metallicity of [Fe/H] = —0.7, 
motivated by results for the disk metallicity distribution presented in 108 and the fact that 
SDSS data are too shallow to include a large fraction of red halo stars. A slightly better 
approach would be to use the disk metallicity distribution from 108 to solve for best-fit 
distance iteratively, but the resulting changes in the photometric distances are negligible 
compared to other systematic errors. 



2.3. The SDSS-POSS Proper Motion Catalog 



Munn et al. 



(I2OO4J ) catalog (distributed 



We take proper-motion measurements from the 
as a part of the public SDSS data releases), which is based on a comparison of astrometric 
measurements between SDSS and a collection of Schmidt photographic surveys. Despite the 
sizable random and systematic astrometric errors in the Schmidt surveys, the combination of 
a long baseline (~ 50 years for the POSS-I survey), and a recalibration of the photographic 
data using the positions of SDSS galaxies (see Munn et al. for details), results in median 
random proper-motion errors (per component) of only ~ 3 mas yr^^ for r < 18 and 
~ 5 mas yr^^ for r < 20. As shown below, systematic errors are typically an order of 
magnitude smaller. At a distance of 1 kpc, a random error of 3 mas yr^^ corresponds 
to a velocity error of ~ 15 km s~^, which is comparable to the radial velocity accuracy 
delivered by the SDSS stellar spectroscopic survey (~ 5.3 km s^^ at g = 18 and 20 km s~^ 
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a.t g = 20.3; Schlaufman et al. 2009). At a distance of 7 kpc, a random error of 3 mas yr~^ 
corresponds to a velocity error of 100 km s^^, which still represents a usable measurement 
for large samples, given that systematic errors are much smaller (~ 20 km s~^ at a distance 
of 7 kpc). The small and well- understood proper-motion errors, together with the large 
distance limit and sample size (proper-motion measurements are available for about 38 
million stars with r < 20 from SDSS Data Release 7) make this catalog an unprecedented 
resource for studying the kinematics of Milky Way stars. 

We warn the reader that proper-motion measurements made publicly available prior to 
SDSS Data Release 7 are known to h ave significant sys tematic errors. Here we use a revised 



set of proper-motion measurements (IMunn et al. 



20081 ). which are publicly available only 



since Data Release 7. As described in the next section, we can assess the error properties of 
this revised proper motion catalog using objects with known zero proper motion - that is, 
distant quasars. 



2.3.1. Determination of Proper Motion Errors Using Quasars 



All known quasars are sufficiently distant that their proper motions are vanishingly 



small compared to the expected random and systematic errors in the Munn et a 



cata log 



20071 ) 



The large number of spectroscopically-confirmed SDSS quasars (iSchneider et al. 
which were not used in the recalibration of POSS astrometry can therefore be used to derive 
robust independent estimates of these errors. In SDSS Data Release 7, there are 69, 916 
quasars with 14.5 < r < 20, redshifts in the range 0.5 < z < 2.5, and available proper 
motions (see Appendix for the SQL query used to select and download the relevant data 
from the SDSS CAS). Within this sample of quasars, the proper motions have a standard 
deviation of ~ 3.1 mas yr~^ for each component (determined from the inter-quartile range), 
with medians differing from zero by less than 0.2 mas yr^^. The dependence of the random 
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error on r-band magnitude is well-described by 

= 2.7 + 2.0 X 100-4 ('■-20) ^as yr"^ (1) 

fitting only to quasars in the range 15 < r < 20. When the measurements of each 
proper- motion component are normalized by a^, the resulting distribution is approximately 
Gaussian, with only ~ 1.8% of the quasar sample deviating by more than 3a from zero 
proper motion. In addition to their dependence on magnitude, the random proper-motion 
errors also depend on position on the sky, but the variation is relatively small (~ 20%, see 
right panels in Figured]). Finally, we find that the correlation between the errors in the two 
components is negligible compared to the total random and systematic errors. 

The median proper motion for the full quasar sample is ~ 0.2 mas yr"^, but the 
systematic errors can be larger by a factor of 2 — 3 in small sky patches, as illustrated in 
Figure [TJ We find that the distribution of systematic proper-motion errors in ~ 100 deg^ 
patches of sky has a width of ~ 0.67 mas yr"^ in each component, about twice as large 
as that expected from purely statistical noise (per bin, using equation ([T])). As the figure 
shows, a few regions of the sky have coherent systematic errors at a level ~ 1 mas yr^^ 
(e.g., the median fii towards / ~ 270°, or fij, towards the inner Galaxy). Therefore, the 
kinematics measured using proper motions in these regions should be treated with caution. 

The largest systematic errors, ~ 1 mas yr"^ for fii, are seen toward / ~ 270° in the 
top left panel in Figure [H which corresponds to 5 < 10°. In this region, the systematic 
deviation of quasar proper motions from zero is approximately parallel to lines of constant 
right ascension, suggesting that the data may be suffering from systematic effects due to 
atmospheric refraction and spectral differences between quasars and galaxies used in the 
recalibration of POSS astrometry. This effect would be strongest for observations obtained 
at high airmass, as are typical for fields at low declination (the POSS data were obtained 
at a latitude of +33°). We find that the median quasar proper motion in the 6 direction is 
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well- described by 

(/i^) = -0.72 + 0.019 S mas yr"^ (2) 
for —5° < S < 30°, where S is in degrees. At 5 > 30", we find (/i^) < 0.2 mas yr~^. 

The observed direction and magnitude of this systematic offset (corresponding to 
an astrometric displacement of up to ~ 30 mas) are consist ent with detailed studies 



200% . 



of atmospheric dispersion effects on observations of quasars ( iKaczmarczik et al. 
Therefore, it is possible that the systematic errors in stellar proper motions (whose spectral 
energy distributions differ less from galaxy spectral energy distributions than is the case for 
quasars) are smaller than implied by Figure [H Nevertheless, we will conservatively adopt 
the quasar proper-motion distributions as independent estimates of systematic and random 
proper motion errors for stars analyzed in this work. In particular, we adopt 0.6 mas yr~^ 
as an estimate for typical systematic proper-motion error. 

The quasar sample has a much narrower color range than that seen in main-sequence 
stars (96% of the quasar sample satisfies —0.2 < g — r < 0.6), and provides a better estimate 
of systematic proper-motion errors for the blue stars than for the red stars. Within the 
above well-sampled color range, we find a median proper-motion gradient with respect to 
the g — r color of < 0.1 mas yr~^ mag~^ (per component). When the fit is extended to 
g — r < 1.6 (using a much smaller number of quasars), the gradient is still smaller than 
0.5 mas yr~^ mag^^. Hence, the proper-motion systematics have a color dependence that is 
smaller than, or at most comparable to, their dependence on sky position. 

A systematic error in proper motion of 0.6 mas yr~^ corresponds to a systematic 
velocity error of 3 km s~^ at 1 kpc, and ~ 20 km s~^ at 7 kpc. In addition, the ~ 5% 
systematic distance errors discussed in § 12.11 are responsible for a ~ 5% systematic 
velocity uncertainty. Hence, for a disk-like heliocentric tangential velocity of 20 km s~^, 
proper-motion systematics dominate at distances beyond ~ 1 kpc. Similarly, for a halo-like 
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heliocentric tangential velocity of 200 km s~^, proper motion systematics will dominate 
at distances greater than 7 kpc. At smaller distances, the dominant systematic in our 
tangential-velocity estimates comes from systematic distance errors. For most of the Galaxy 
volume analyzed in this work, the systematic distance errors dominate over systematic 
proper-motion errors. 



2.4. Comparison of Proper Motions with Independent Measurements 



As further tests of the proper-motion errors, we have analyzed two independent sets 
of measurements. As shown below, they confirm the results based on our analysis of the 
quasar sample. 

We have compa red the SDSS-POSS proper motions to proper-motion measurements 



by iMajewskil (jl992l ) for a sample of 326 stars observed towards the North Galactic 
Pole. The measurements in the Majewski sample have random errors that are three 
times smaller, and comparable, but most likely, different systematic errors. The median 
proper- motion differences between the two data sets are below 1 mas yr~^, with an rms 



scatter 3 — 4 mas yr ^ (per coordinate). Hence, this comparison is cons istent wit 



h our error 



estimates discussed in the preceding section, and with the estimates of 



Maiewskil fll992h 



We have also compar ed the SDSS-POSS pr oper motions with proper motions from the 



SDSS stripe 82 region. In 



Bramich et al. 



(120081 ) ■ proper motions are computed using only 



SDSS data, and thus they are expected to have different, and probably smaller, systematic 
errors than the SDSS-POSS proper motions (random errors for the stripe 82 proper motions 
are larger by about a factor of two). For ~500,000 stars with both SDSS-POSS and 
Bramich et al. proper-motion measurements, we find the median differences and the rms 
scatter to agree with expectation. A single worrisome result is that the median difference 



- 15 - 



between the two data sets is a function of magnitude: we find a gradient of 0.8 mas yr~^ 
between r = 15 and r = 20. It is more likely that this gradient is due to systematic errors 
in centroiding sources on photographic plates, rather than a problem with SDSS data. 
This gradient corresponds to a systematic velocity error as a function of distance, Av ~ 4 
(D/kpc) km s~^. For example, a halo star at 5 kpc, with a relative velocity of 200 km s~^, 
would have a systematic velocity uncertainty of 10%. This systematic error is comparable 
to other sources of systematic errors discussed above, and has to be taken into account 
when interpreting our results below. 

2.5. The Main Stellar Samples 

When using proper motions, random errors in the inferred velocities have a strong 
dependence on magnitude, and therefore distance, while systematic errors are a function 
of position on the sky, as discussed above. Random errors in radial-velocity measurements 
also depend on magnitude, as fits to spectral features become more difficult at lower 
signal-to-noise ratios. As such, when radial-velocity and proper-motion measurements are 
analyzed simultaneously, the systematic and random errors combine in a complex way - 
care is needed when interpreting the results of such an analysis. 

In order to minimize these difficulties, we separately analyze the proper-motion sample 
and the much smaller sample of stars with radial velocities. Furthermore, motivated by the 
metallicity distribution functions presented in 108, we separately treat the low-metallicity 
"halo" stars and the high-metallicity "disk" stars. For these two samples, wc require 
g — r < 0.6, the regime in which the photometric metallicity estimator is believed to be 
accurate. Finally, we discuss a sample of "red" stars with — r > 0.6 (roughly, g — i > 0.8), 
which are dominated by nearby (< 2 kpc) disk stars. 
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These samples are selected from SDSS Data Release 7 using the following common 
criteria: 

1. Unique unresolved sources that show subarcsecond parallax: binary processing 
flags DEBLENDED_AS_MOVING, SATURATED, BLENDED, BRIGHT, and 
NODEBLEND must be false, and parameter nCHILD=0 

2. The interstellar extinction in the r band. A,. < 0.3 mag 

3. Dust-corrected magnitudes in the range 14.5 < r < 20 mag 

4. High galactic latitudes: \b\ > 20° 

5. Proper motion available. 



yieldin g 20.1 million stars. The dust corrections, Ar, were computed using the lSchlegel et al. 



(119981 ) dust maps, with conversion coefficients derived assuming an Ry = 3.1 dust model. 



The intersection of the following color criteria then selects stars from the main stellar locus: 

• Blue stars (6.9 million): 

1. 0.2 <{g-r) < 0.6 

2. 0.7 <{u-g) < 2.0 and -0.25 < [g - r) - 0.5(u - g) < 0.05 

3. -0.2 < 0.35(^ -r)-{r-i) < 0.10 

• Red stars (11.9 million): 

1. 0.6 <{g-r) < 1.6 

2. -0.15 < -0.270r + O.SOOi - 0.534z + 0.054 < 0.15, 
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where the last condition is based on a "principal colo r" orthogon a. 



Ivezic et al 



to the stellar locus in 



fl2004f ). This condition 



the i — z Ys. r — i color-color diagram, as defined in 
allows for a 0.15 mag offset from the stellar locus. 

During the analysis, "blue" stars are often further divided by photometric metallicity 
(see below for details) into candidate "halo" stars ([Fe/H] < —1.1) and candidate "disk" 
stars ([Fe/H] > —0.9). Subsamples with intermediate metallicities include non-negligible 
fractions of both halo and disk stars. Although the reduced proper-motion diagram 
is frequently used for the separation and analysis of these two populations, we find it 
inadequate for our purposes; the vertical gradient in rotational velocity blurs the kinematic 
distinction between disk and halo (for a discussion, see SIJ08), and thus this method is 
applicable only to stars with significant proper motion (leading to severe selection effects). 
Although metallicities are not available for red stars, results from 108 imply that they are 
dominated by the disk population (red stars can only be seen out to ~ 2 kpc). 

For each of the subsamples defined above, we further separate those objects with SDSS 
spectroscopic data (see Appendix for a sample SQL query) into independent subsamples. 
In total, these spectroscopic subsamples include 172, 000 stars (out of 352, 000 stars with 
spectra), after an additional requirement to select only main-sequence stars; that is, stars 
with \og{g) > 3.50. Of the stars with spectroscopic data. 111, 000 are blue (0.2 < g — r < 0.6) 
and 61,000 are red (0.6 < g — r < 1.6). When separati ng low- and high-me t allicity stars 



AUende Prieto et al 



(120061) for 



with spectra, we use the spectroscopic metallicity (see 
details). Due to increased difficulties with measuring [Fe/H] for red stars {g — r > 0.6) 
from SDSS spectra, we adopt [Fe/H] = —0.7 for all such stars; this value is the median 
spectroscopic [Fe/H] for stars with 0.6 < 5^ — r < 1.3 (a=0.4 dex). For over 90% of ~30,000 



^Note that the majority of stars with g — r > 1.2 do not have reliable measurements of 
log{g) - we assume that all stars with g — r > 1.2 are main-sequence stars. 
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stars with g — r > 1.3, [Fe/H] is not successfully measured. 



2.6. Coordinate Systems and Transformations 

Following J08 and 108, we use a right-handed, Cartesian Galactocentric coordinate 
system defined by the following set of coordinate transformations: 

X^Rq-D cos(Z) cos(6) 

Y^-D sin(/) cos(6) (3) 
Z sin(6), 

where Rq = 8 kpc is the adopted distance to the Galactic center, D is distance of the 
star from the Sun, and (Z, b) are the Galactic coordinates. Note that the Z — plane 
passes through the Sun, not the Galactic center (see J08), the X axis is oriented toward 
/ = 180°, and the Y axis is oriented toward / = 270° (the disk rotates toward / ~ 90°). 
The main reason for adopting a Galactocentric coordinate system, rather than a traditional 
heliocentric system, is that new data sets extend far beyond the solar neighborhood. 

We also employ a cylindrical coordinate system defined by 

(4) 



R = VX^ + Y\ 4> = tan-^ (J^^ . 



The tangential velocity, v, is obtained from the proper motion, /x, and the distance D 

by 

V = 4.74 — - km s-^ (5) 

mas yr~^ kpc 

Given the line-of-sight radial velocity, I'^ad) and the two components of tangential 
velocity aligned with the Galactic coordinate system, vi and v^, the observed heliocentric 
Cartesian velocity components are given by: 

■^x ^ — ~Vrad cos(Z) cos(&) + Vb cos(Z) sin(6) + vi sin(Z) 
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obs 

-'y 



,obs 



-i^rad sin(/) cos(6) + Vh sin(/) sin(6) — vi cos(/) 
-t;radsin(6) + VbCos{b). 



(6) 



These components are relat ed to the traditional UVW nomenclature by, vx 



Vy 



-V, and vz = W, e.g., 



Binnev fc Merrifieldl (119981 ). 



In order to obtain the Galactocentric cylindrical velocity components, we must first 
correct for the solar motion. Ta king into account HI measurements of the Galactic 
rotation curve (( 



Gunn et al 



(IFeast fc Whitelock 



19791 ) and Hipparcos measurements of Cepheid proper motions 



19971 ), we adopt wlsr = 220 km s ^ for the motion of the local stand ard 



of rest and R, 



■0 



kpc (for an analysis of other recent measurements, see 



Bovy et al. 



20091 ). The adopted value of Rq is motivated by geometrical 



of stars around Sgr A*, which yield Rq = 7.94 ± 0.42 kpc (lEisenhauer et al. 



measurements o: 



the motions 



20031 ). For 



the solar peculiar motion, we adopt the HIPPARCOS result, Vy 



0,pec 



-10.0 ±0.4 km s-\ 



0,pec _ 

Uy — 

also see 



—5.3 ± 0.6 km s ^, and 



Hogg et al. 



0,pcc 



7.2 ± 0.4 km s^^ flDehnen fc BinnevI (119981 ) 



(120051 )). Using these values, along with equation ([6]), we obtain the 



Galactocentric velocity components: 



(7) 



where 

,0 



-10 km s ^, Vy 



-225 km s ^, and ff = 7 km s (note that 



Vy = -VlSR + Vy 



0,pec\ 



Below, we discuss attempts to directly determine the solar peculiar 



motion (§ 14.21) and f lsr (§ 15.31) from our data. 



Finally, the cylindrical components, vji and v^, can be computed using a simple 
coordinate rotation, 

X Y 

vr = vx-p: + Vy — 
ri ri 

Y X 
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Note that, in our adopted system, the disk has a prograde rotation = —220 km s~^; 
retrograde rotation is indicated by > 0. Stars with vr> move away from the Galactic 
center, and stars with vz > move toward the North Galactic Pole. 



2.7. Analysis Philosophy 

Such a massive data set, extending to a large distance limit and probing a large fraction 
of the Galaxy volume, can be used to map the kinematics of stars in great detail. It can also 
be used to obtain best-fit parameters of an appropriate kinematic model. However, it is not 
obvious what model (functional form) to chose without at least some preliminary analysis. 
Hence, in the next two sections, we first discuss various projections of the mult i- dimensional 
space of the available observable quantities and obtain a number of constraints on the 
spatial variation of stellar kinematics. We then synthesize all of the constraints into a model 
described in § Before proceeding with our analysis, we provide a brief summary of the 
first two papers in this series, whose results inform our subsequent analysis. 



2.8. A Summary of Papers I and II 



Using photometric data fo r 50 million stars from SDSS Data Release 4 



( lAdelman-McCarthy et al.ll2006l ). sampled over a distance range from 100 pc to 15 
kpc, JOS showed that the stellar number density distribution, p{R, Z, 0) can be well- 
described (apart from local overdensities; the JOS best-fit was refined using residual 
minimization algorithms) as a sum of two cylindrically symmetric components: 



p{R, Z, 



Pn{R, Z)+ph{R, Z). 



(9) 
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The disk component can be modeled as a sum of two exponential disks 

/ \Z+Zq\ (R-Rq) \Z+Zq\ (R-Rq) \ 

PD{R,Z)^PD{Re){c «i ^1 +eDe «2 ^2 , (10) 



while the halo component requires an oblate power-law model 

The best- fit parameters are discussed in detail by JOS. We have adopted the following 
values for the parameters relevant to this work (second column in Table 10 from JOS): 
^0=25 pc, Hi = 245 pc, H2 = 743 pc, en = 0.13, €h = 0.0051, qh = 0.64, and 
uh = 2.77. The normalization pd{Rq) (which is essentially the local luminosity function 
for main-sequence stars) is listed in JOS as a function of color. 

Using a photometric metallicity estimator for F/G main-sequence stars, 108 obtained 
an unbiased, three-dimensional metallicity distribution of ~2.5 million F/G stars at 
hehocentric distances of up to ~S kpc. They found that the metallicity distribution 
functions (MDF) of the halo and disk stars are clearly distinct. The median metallicity of 
the disk exhibits a vertical (with respect to the Galactic plane, Z) gradient, and no gradient 
in the radial direction (for Z > 0.5 kpc and 6 < < 10 kpc). 

Similarly to the stellar number density distribution, p{R,Z), the overall behavior of 
the MDF p([Fe/H]|it!, Z) for disk stars can be well-described as a sum of two components 

p{x = [Fe/R]\R, Z, 0) = [1 - MR, Z)] pd{x\Z) + MR, Z) p„{x), (12) 

where the halo star-count ratio is simply, 

■""^ ' ' pd{R,Z)+ph{R,Z) ^ ' 

The halo metallicity distribution, p^([Fe/H]), is spatially invariant within the probed 
volume, and well-described by a Gaussian distribution centered on [Fe/H] = —1.46, with an 
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intrinsic (corrected for measurement errors) width an — 0.30 dex. For \Z\ < 10 kpc, an 
upper limit on the halo radial metallicity gradient is 0.005 dex/kpc. 

The disk metallicity distribution varies with Z such that its shape remains fixed, while 
its median, /xd, varies as 

/^d(^) =/^oo + A^e-^, (14) 

with the best-fit parameter values i?^ = 0.5 kpc, //oo = —0.82, and = 0.55. The shape of 
the disk metallicity distribution can be modeled as 

Pd{x^ [Fe/H]|Z) = 0.63G'[x|// = a(Z),(7 = 0.2] + 

Q.37G[x\n = a{Z) + O.U,a = 0.2], (15) 

where the position a and the median are related via a{Z) = fJ^oiZ) — 0.067 (unless 
measurement errors are very large). 

The main result of this third paper in the series is the extension of these results for 
number density and metallicity distributions to include kinematic quantities. 



3. ANALYSIS OF THE PROPER MOTION SAMPLE 

We begin by analyzing the proper-motion measurements of stars observed toward the 
North Galactic Pole. In this region, the Galactocentric azimuthal velocity, v,f,, and radial 
velocity, vr, can be determined directly from the proper-motion measurements (that is, 
without knowledge of the spectroscopically-determined radial velocity, Vrad)- In this way, 
we can study the kinematic behavior of stars as a function of metallicity and distance from 
the Galactic plane, Z. We then extend our analysis to the entire meridional Y = plane, 
and study the variation of stellar kinematics with R and Z. In the following section, we 
only consider the northern Galactic hemisphere, where most of the proper-motion data are 
available. 
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3.1. Kinematics Towards the North Galactic Pole 

We select three stellar subsamples in the region b > 80°, including 14, 000 blue disk 
stars at Z < 7 kpc, 23, 000 blue halo stars at Z < 7 kpc, and a sample of 105, 000 red stars 
at Z < 1 kpc. In Figure [2], we plot the distribution of vs. vr for ~ 6, 000 blue disk and 
halo stars at Z = 4 — 5 kpc. In this and all subsequent two-dimensional projections of the 



velocity distribution we plot smoothed, color-coded maps, where t 



are estimated using the Bayesian density estimator of llvezic et al. 



le ve locity distributions 



(12005 



see their Appendix 



for the derivation and a discussion). At a given position, the density is evaluated as 

where di is the distance in the velocity-velocity plane, and we sum over the = 10 nearest 
neighbors. The normalization constant, C, is easily evaluated by requiring that the density 
summed over all pixels is equal to the total number of data points divided by the total area. 
The grid size is arbitrary, but the map resolution is determined by the density of points - 



we ch oose pixel size equal to one half of the mean velocity error. As shown by llvezic et al. 



( I2OO5I ). this method is superior to simple Gaussian smoothing. For comparison, we also plot 



linearly-spaced density contours. 

The six panels of Figure [2] demonstrate the variation of kinematics with metallicity, 
with the full range of metallicities (—3 < [Fe/H] < 0) plotted in the upper left panel and 
subsamples with increasing metallicity running from left-to-right, top-to-bottom. The mean 
azimuthal velocity varies strongly with metallicity, from a non-rotating low-metallicity 
subsample with large velocity dispersion (top center panel) to a rotating high-metallicity 



sample with much smaller dispersion (bottom right panel). This strong meta 
correlation is qualitatively the same as discussed in the seminal paper by 



licity-kinema tic 



Eggen et al. 



( 19621 ) . except that here it is reproduced in situ with a ~100 times larger, nearly-complete 



sample, thus extending it beyond the solar neighborhood. There are some indications of 
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substructure in the velocity distribution, but much of it remains unresolved due to the large 
velocity measurement errors. 

The substructure becomes more apparent in Figure [3], where we plot the same 
velocity-space projection for 60, 000 stars within Z < 2.5 kpc. In this figure, the panels show 
subsamples of increasing distance from the Galactic plane, beginning with Z = 0.1 — 0.2 kpc 
in the upper left panel (note the changing axes between the top and bottom rows). 
The substructure seen in the closest bin probed by red stars is v ery similar to the 



19981 ). These results were 



substructure seen in the local HIPPARCOS sample (iDehnenl 
based on a maximum-likelihood analysis over the entire sky, while our result arises from 
a direct mapping of the velocity distribution of stars selected from a small region (~ 300 
deg^). Using a subsa mple of ~ 17, OOP HIPP ARCOS stars with full three-dimensional 



veloci t y info rmat ion 



Nordstrom et al 



(12007 



(120041 ): 



Famaey et al. 



(!2005h : 



Holmberg et al. 



20091 ) have detected the same kinematic morphology. The similarity between these 



HIPPARCOS-based velocity d istributions and ours, including the multi-modal behavior 



reminiscent of moving groups ( lEggen 



19961 ). is quite encouraging, given the vastly different 
data sources. The similarity of observed substructure with moving groups is even more 
striking for stars from a closer distance bin (Z = 50 — 100 pc ), matched to distances probed 



by the HIPPARCOS sample (see Figure Hj). As suggested by 



De Simone et al. 



(120041 ) ■ these 



moving groups may arise from irregularities in the Galactic gravitational potential. 

The remainder of our analysis will focus on blue stars, which sample a much larger 
distance range. For a detailed study of the velocity distribution of nearby red stars, 
including a discussion of non-Gaussianity, vertex deviations, and difficulties with traditional 
thin/thick disk separation, we refer the reader to Kowalski et al. (in prep.). 

The dependence of the rotational velocity on height above the Galactic plane is shown 
in Figure [51 The two subsamples display remarkably different kinematic behavior (first 
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seen locally by lEggen et al.l Il962l ) with halo stars exhibiting a small constant rotational 
motion (~ —20 km s^^), and disk stars exhibiting a large rotational- velocity component 
(~ —200 km s~^ at Z ~ 1 kpc) that decreases with height above the Galactic plane. 

We have performed the same analysis using proper motions based only on POSS data, 
with SDSS positions not included in the proper-motion fit (not publicly availablcl). While 
random proper-motion errors become larger when SDSS data are not used, the median 
rotational velocity for halo stars decreases to only 5 km s~^, suggesting that the apparent 
rotational motion in the halo subsample is influenced by systematic errors. These tests 
also suggest that the leading contribution to systematic proper-motion errors could be 
a difference between the SDSS (digital data) a nd POSS (dig i tized photographic data) 



centroid-determination algorithms. In addition. 



Smith et al. 



(120091) did not detect halo 



rotation using a smaller s ample, but with more robus t proper motion measurements based 



on only SDSS data, while 



AUende Prieto et al 



( 120061 ) found no evidence for halo rotation 



using SDSS DR3 radial velocities. We conclude that the net halo rotation in the direction 
of the North Galactic Pole is \v^ot\ ^ 10 km s^^. In addition, the measured halo velocity 
dispersion increases with Z, but when random measurement errors are taken into account, 
the data are consistent with a constant dispersion of a^j^ = 85 ± 5 km s^^ (derived using the 
test described in § |5]) . 

The decrease of rotational velocity with Z for disk stars (often referred to as asymmetric 
drift, velocity lag, or velocity shear; see §3.4 of 108 for more details and references to related 
work) is in agreement with a preliminary analysis presented in 108. We find that the 
observed behavior in the Z = 1 — 4 kpc range can be described by 

ry 1.25 



-205 + 19.2 



kpc 



km s ^. 



;i7) 



^Available from J. Munn on request 
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The measured rotational velocity dispersion of disk stars increases with Z faster than can 
be attributed to measurement errors. Using a functional form a = a + h\ZY^ we obtain an 
intrinsic velocity dispersion fit of 

ry 2.0 



= 30 + 3.0 



kpc 



This function and the best-fit rotational velocity for halo stars are shown as dotted 
lines in the bottom right panel of Figure [5] (see Table 1 for a summary of all best-fit 
parameters). 108 fit a linear model to vs. Z, but the difference between this result and 
their equation (15) never exceeds 5 km s~^ for Z < 3 kpc. The errors on the power-law 
exponents of equations (fTTI) and (fT8|) are ~ 0.1 and ~ 0.2, respectively. 

However, a description of the velocity distribution based solely on the first and second 
moments (equations (|T71) and (|T8l) ) does not fully capture the detailed behavior of our data. 
As already discussed by 108, the rotational- velocity distribution for disk stars is strongly 
non-Gaussian (see their figure 13). It can be formally described by a sum of two Gaussians, 
with a fixed normalization ratio and a fixed offset of their mean values for |Z| < 5 kpc, 

Pd{x = v^\Z) = 0.75 G[x\vrXZ),ai] + 0.25 G[x\vn{Z) - 34km s-\a2\, (19) 

where 



Vn{Z) = -194 + 19.2 



kpc 



1.25 



km s-\ (20) 



The intrinsic velocity dispersions, ai and (J2, are modeled as a + b\Z\'^, with best-fit 
parameters listed in Tabled] (see and cr^). Closer to the plane, in the 0.1 < Z < 2 kpc 
range probed by red stars, the median rotational velocity and velocity dispersion are 
consistent with the extrapolation of fits derived here using much more luminous blue stars. 

Figure |6] shows the f,^ distribution for four bins in Z (analogous to figure 13 from 
108), overlaying two-component Gaussian fits with the measurement errors and Vn{Z) as 
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free parameters. The mean velocity and velocity dispersion exhibit ~ 10 km s~^ variations 
relative to their expected values; while such deviations could be evidence of kinematic 
substructure, they are also consistent with the plausible systematic errors. We conclude 
that equations (fT9|) and (!20|) provide a good description of the disk kinematics for stars 
observed toward the North Galactic Pole, within the limitations set by the random and 
systematic errors in our data set. 

The fits to the observed velocity distributions for halo and disk stars are shown in 
Figure M and demonstrate that the vertical gradients in median rotational velocity and 
velocity dispersion for disk stars seen in Figure are not due to contamination by halo 
stars. To quantitatively assess the impact of "population mixing" as a result of the adopted 
metallicity-based classification ([Fe/H] < —1.1 for "halo" stars and [Fe/H] > —0.9 for 
"disk" stars) on our measurements of these gradients, we have performed a series of Monte 
Carlo simulations. Assuming that the fits shown in Figure |6] accurately depict the intrinsic 
velocity distributions, and adopting analogous fits for their metallicity distributions (108, 
their Figure 7), we have estimated the expected bias in median velocity and velocity 
dispersion for each population as a function of their relative normalization. As shown in 
Figure 6 from 108, the fraction of halo stars increases with distance from the plane, from 
about 0.1 at Z = 1 kpc to about 0.9 at Z = 5 kpc. We find that the median rotation 
velocity and velocity dispersion biases are < 10 km s^^ for disk stars at Z < 3 kpc, as well 
as for halo stars at Z > 3.5 kpc. Furthermore, the biases are < 20 km s~^ for Z < 4 kpc 
for disk stars and to Z > 3 kpc for halo stars. At Z = 3 kpc, the contamination of 
both disk and halo subsamples by the other population is typically ~ 10 — 15%. This 
small contamination and the use of the median (as opposed to the mean) and dispersion 
computed from the interquartile range results is reasonably small biases. With the adopted 
metallicity cutoffs, the sample contamination reaches 50% at Z = 1 — 1.5 kpc for halo stars 
and at Z = 4.5 kpc for disk stars. 
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The dependence of the Galactocentric radial velocity on Z is shown for halo and disk 
subsamples in Figure [71 The median values (bottom left panel) are consistent with zero, 
within the plausible systematic errors (10 — 20 km s""^), at all Z . The intrinsic dispersion 
for halo stars is consistent with a constant value of = 135 ± 5 km s~^. For disk stars, 
the best-fit functional form a = a + b\Z\'^ is 

1.5 

km s-\ (21) 

The (T^/cr^ ratio has a constant value of ~ 1.35 for Z < 1.5 kpc, and decreases steadily at 
larger Z to ~ 1 at Z ~ 4 kpc. 



= 40 + 5 



'R 



kpc 



3.2. Kinematics in the Meridional F ~ Plane 

The analysis of the rotational-velocity component can be extended to the meridional 
plane defined by F = 0, for which the longitudinal proper motion depends only on the 
rotational- velocity component and the latitudinal proper motion, f;,, is a linear combination 
of radial and vertical components, 

Vh = sm{b)vR + cos{b)vz- (22) 

Figure [8] plots and Vb as functions of R and Z for halo and disk stars within 10° of the 
meridional plane. The median Vb is close to zero throughout most of the plotted region, as 
would be expected if the median vr and vz are zero (the behavior of vz is discussed in the 
next section). One exception is a narrow feature with Vb ~ —100 km s^^ for i? < 4 kpc. 
While a cold stellar stream could produce such a signature, its narrow geometry points 
directly at the observer. This behavior is also consistent with a localized systematic 
proper-motion error. Indeed, the bottom left panel in Figure [T] shows that the systematic 
latitudinal proper-motion error at / ~ 0°, 6 ~ 45° is about 1 mas yr~^, corresponding to a 
velocity error of ~ 100 km s~^ at a distance of 7 kpc. 
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As seen in the upper left panel of Figure [H the median for halo stars is close to 
zero for R < 12 kpc. In the region with R> 12 kpc and Z < 6 kpc, the median indicates 
a surprising prograde rotation in excess of 100 km s^^. This behavior is also seen in disk 
stars, and is likely due to the Monoceros stream, which has a metallicity intermediate 
between disk and halo stars and rotates faster than disk stars (see §3.5.1 and §3.5.2 in 108). 
There is also an indication of localized retrograde rotation for halo stars with Z ~ 9 kpc 
and ~ 15 kpc (corresponding to I ~ 180°, h ~ 50°, and a distance of ~ 11 kpc). Stars 
with Z = 8 — 10 kpc and i? = 15 — 17 kpc have median larger by 40 km s~^ (a ~ la 
effect) and median [Fe/H] larger by 0.1 dex (~ 5cr effect) than stars with Z = 8 — 10 kpc 
and i? = 7 — 13 kpc. A systematic error in /i^ of ~ 0.8 mas yr~^ is required to explain 
this kinematic feature as a data problem (although this would not explain the metallicity 
offset). However, the top right panel in Figure [T] shows that the systematic /i; errors in this 
sky region are below 0.5 mas yr^^, so this feature may well be real. We note that in roughly 



Grillmair &: Dionatos 



(120061) have 



the same sky region and at roughly the same distance, 
detected a narrow stellar stream. 

In order to visualize the extent of "contamination" by the Monoceros stream, we 
replace the rotational velocity for each disk star by a simulated value drawn from the 
distribution described by equation ( IT9l) . We subtract this model from the data, and 
the residuals are shown in the right-hand panel of Figure |9l The position of the largest 
deviation is in excellent agreement with the position of Monoceros stream quantified in 108 
{R = 15 — 16 kpc and \Z\ ~ 3 — 5 kpc). Further evidence for the presence of the Monoceros 
stream is shown in Figure [101 in which we analyze v^f, vs. [Fe/H], as a function of i?, for blue 
stars at Z = 4 — 6 kpc. As is evident in the bottom right panel, there is a significant excess 
of stars at i? > 17 kpc with —1.5 < [Fe/H] < —0.5 that rotate in a prograde direction with 
~ 200 km s-^ 
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4. ANALYSIS OF THE SPECTROSCOPIC SAMPLE 

Despite its smaller size, the SDSS DR-7 spectroscopic sample of ~ 100, 000 main- 
sequence stars is invaluable, because it enables a direcio study of the three-dimensional 
velocity distribution. The sample extends to a distance of ~ 10 kpc, at which it can 
deliver velocity errors as small ~ 10 km s~^ (corresponding tangential velocity errors are 
~ 150 km s~^ at a distance of 10 kpc). For each object in the SDSS spectroscopic survey, 
its spectral type, radial velocity, and radial-velocity error are determined by matching the 



measured spec trum to a set o 



stellar library ( iPrugniel et al 



stel 



ar templates, which were calibrated using the ELODIE 
20071 ). Random errors on the radial- velocity measurements 
are a strong function of spectral type and signal-to-noise ratio, but are usually < 5 km s~^ 
for stars brighter than g ~ 18, rising sharply to ~ 15 km s~^ for stars with r = 20. We 
model the behavior of the radial-velocity errors as 



a,ad = 3 + 12x 10°-^(^'-20) kjns-i. 



(23) 



We begin our analysis with blue disk and halo stars, and then briefly discuss the 
kinematics of nearby red M stars. 



4.1. Velocity Distributions 

We select 111,000 stars with 0.2 < ^ - r < 0.6 (74,000 have b > 20°) and, using 
their spectroscopic metallicity, separate them into 47, 000 candidate halo stars with 
[Fe/H] < —1.1, and ~ 53,000 disk stars with [Fe/H] > —0.9. Assuming the spectroscopic 



■^Statistical de proiection rn e thods , such as that recently applied to a subsample of M 
stars discussed by iFuchs et al.l (120091 ). can be used to indirectly infer the three-dimensional 
kinematics from proper motion data. 
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metallicities accurately separate disk from halo stars, the use of photometric metallicity for 
the same selection would result in a contamination rate of 6% for the halo subsample, and 
12% for the disk subsample. 

The dependence of the median vertical velocity, Vz-, and its dispersion on height above 
the disk, is shown in Figure [TT] for the halo and disk subsamples. The median values 
of Vz are consistent with zero to better than 10 km s~^ at Z < 5 kpc, where statistical 
fluctuations are small. 

As with and a^, the vertical velocity dispersion can be modeled using a constant 
dispersion for halo stars (a^ = 85 km s^^), while for disk stars, the best-fit functional form 
is 



= 25 + 4 



z 



km s-^ (24) 



kpc 

The other two velocity components behave in a manner consistent with equations ([1] 
( !T8|) . and (!2T|) . just as they did in the proper- motion sample. This is encouraging, because 
the spectroscopic sample is collected over the entire northern hemisphere, unlike the 



proper-motion subsample studied in § 13. H which is limited to 6 > 80°. 

The availability of all three velocity components in the spectroscopic sample makes 
it possible to study the orientation of the halo velocity ellipsoid. Figure [T2] shows 
two-dimensional projections of the velocity distribution for subsamples of candidate halo 
stars with 0.2 < g — r < 0.4. The top row corresponds to stars above the Galactic plane 
at 3 < Z/kpc < 4, while the bottom row is for stars the same distance below the plane. 
The velocity ellipsoid is clearly tilted in the top- and bottom-left panels, with a tilt angle 
consistent with tan^^{vz/vfi) = R/z. While the tilt-angle errors are too large to obtain 
an improvement over existing measurements of Rq, it is remarkable that the northern and 
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southern subsamples agree so welo In addition, when the Z = 3 — 5 kpc sample is divided 
into three subsamples with 7 < R/kpc < 11, the tilt angle varies by the expected ~ 8° in 
the correct direction (see Figure [T51) . For all bins in the R — Z plane, the best-fit tilt angle 
is statistically consistent (within 5°) with ta.n~^{vz/vii) = R/z. The other two projections 
of the velocity distribution for halo stars do not exhibit significant tilts to within ~ 3° . 

If we transform the velocities to a spherical coordinate system, 

R Z 

Vr = Vr— + Vz — 

ve = - vz^, (25) 

where r = R^^ = {R^ + Z'^Y^'^ is the spherical Galactocentric radius, we find no statistically 
significant tilt in any of the two-dimensional velocity-space projections for halo stars (with 
tilt-angle errors ranging from ~ 1° to ~ 5°). 

As shown in Figure [HI we see no evidence for a velocity ellipsoid tilt in vz vs. vr for 
the disk stars. The plotted subsamples are again selected to have colors 0.2 < g — r < 0.4, 
but are selected closer to the Galactic plane, \Z\ = 1.5 — 2.5 kpc, in order to improve 
statistics and reduce contamination from halo stars. The velocity-ellipsoid tilt is consistent 
with zero within ~ Icr, and alignment of the velocity ellipsoid with the spherical coordinate 
system of equation f l2^ is ruled out at a ~ 2o" or greater confidence level for each of five 
analyzed R — Z bins (i? = 6 — 11 kpc, with AR = 1 kpc). We conclude that there is no 
evidence for a velocity-ellipsoid tilt in the disk subsample, but caution that, due to the small 
Z range, the data cannot easily distinguish between cylindrical and spherical alignment. A 
model prediction for velocity ellipsoid tilt is discussed below. 



^A plausible, if somewhat optimistic, tilt-angle uncertainty of 1° corresponds to an Rq 
)r of 0.5 kpc; extending the sample to \Z\ = 8 kpc could deliver errors of 0.3 kpc per bin 



of a similar size. 
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The vertex deviation is analogous to the velocity-elhpsoid tih discussed above, but 
is defined in the vs. vr plane instead of the vz vs. vr plane. The same plots for red 
{g — r > 0.6, median 1.2) disk stars are shown in the center top and bottom panels in 
Figure [T5l These stars can be traced closer to the plane, \Z\ = 0.6 — 0.8 kpc; in both 
hemispheres, the data are consistent with a vertex deviation of ~ 20°, with an unce r tainty 



of ~ 10°. This result is consistent with the vertex deviation obtained by 



Fuchs et al. 



(120091). 



Another interpretation for the vs. vr distribution of disk stars invokes a two- 
component velocity distribution, which can result in a similar deviation even if each 
component is perfectly symmetric in the cylindrical coordinate system. Kowalski et al. 
(in prep.) find that the vs. vr distribution for red stars toward the North Galactic 
Pole, with 0.1 < Z/kpc < 1.5, can be fit by a sum of two Gaussian distributions that are 
offset from each other by ~ 10 km s~^ in each direction. This offset results in a non-zero 
vertex deviation if the sample is not large enough, or if measurements are not sufficiently 
accurate, to resolve the two Gaussian components. This double- Gaussian structure would 
be at odds with the classical description based on the Schwarzschild approximation - we 
refer the interested reader to the Kowalski et al. study for more details. Unfortunately, the 
spectroscopic samples are not large enough to distinguish a two-component model from the 
standard interpretation. 



4- .1.1. A Model Prediction for Velocity Ellipsoid Tilt 



The tilt of the velocity el 
developed in 



ipsoid tilt in the vz-vr plane is calculated using methods 



Kent fc de Zeeuwl (1l99ll ). who show that the tilt of the ellipsoid at any point 
in the Milky Way depends not only on the gravitational potential of the Galaxy but 
also on the isolating integrals of a particular orbit and thus (weakly) on the distribution 



function of stellar velocities. We use the model of the Milky Way gravitational field from 
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Carlberg &: InnanenI (119871 ). scaled to a solar radius of 8 kpc and a circular velocity of 
220 km s^^. At a particular point in the Galaxy, we integrate eight orbits, each one 
launched from that point, with velocities in each of the R, Z, and (p coordinates that are 
±1(7 from the systemic velocity. We use the observed velocity distributions to obtain the 
systemic velocity and the la values. Furthermore, using the "least-squares fitting" method, 
we determine the parameters of the prolate spheroidal coordinate system that best matches 
each orbit - from this, we can compute its contribution to the velocity ellipsoid. For some 
orbits (especially those of halo stars), the orbits are not local to the solar neighborhood, and 
the local fitting method cannot fit the entire orbit. For these cases, we confine the orbits 
to radii greater than 4 kpc. Since the only purpose of the fit is to model the kinematics in 
a small volume of space, the fact that the fit is no longer global is of no consequence. The 
scatter in tilt angle among the individual orbits at each position depends only weakly on 
position (~ 10%). 

Table [2] gives the predicted tilt angles for each position plotted in Figures [T2lfT5l For 
the halo stars plotted in Figures fT2lfT3l where the ellipsoid shows the most distinctive tilt, 
the predictions match the observed tilts quite well, but are slightly steeper at Z > 3 kpc 
than the Z/R relation predicted for spherical symmetry. For the disk stars in Figures [T^lfT5| 
the predicted tilts lie between the cases of cylindrical and spherical symmetry. The observed 
ellipsoid is sufficiently round, however, that no definitive comparison can be made. 



4.2. Direct Determination of the Solar Peculiar Motion 



If there is no net streaming motion in the Z direction in the solar neighborhood, the 
median heliocentric, v'^^, for nearb y stars should be equal t o ff (7 km s~^, based on an 



analysis of HIPPARCOS results by 



Dehnen fc Binney 



19981 ). We do not expect a large 



velocity gradient within ~ 1 kpc from the Sun, so we select from the spectroscopic sample 
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~ 13, 000 M dwarfs with 2.3 < g — i < 2.8. In the northern hemisphere, we have 5, 700 stars 
with a median hehocentric velocity, {v°^^) = —1.8 km s^^, while for stars in the southern 
hemisphere we obtain {v'^^) = —11.0 km s~^. This difference is likely due to a systematic 
radial-velocity error. If we simultaneously vary an assumed radial-velocity error, Arad, and 
the solar peculiar motion, f^, while requiring that the median f^^*^ should be the same for 
both hemispheres, we obtain Arad = 5.0 ± 0.4 km s^^ and ff = 6.5 ± 0.4 km s~^. This 
value for is in excelle nt agreement with the HIPPARCOS value of 7.2 ± 0.4 km s~^ 



(jPehnen fc Binney 



19981 ). This systematic offset in SDSS radial velocities is probably due 



to the small number of EL ODIE templates 



of improved templates from 



Bochanski et al. 



or red stars. It is likely that the adoption 



(j2007al ) will yield smaller systematic errors. 



A similar analysis for blue stars does not yield a robust detection of the velocity offset. 
A detailed comparison of SDSS radial velocities with radial-velocity standards from the 
literature arrived at the same null result for blue stars. 



of 



As with the vertical component of the solar peculiar motion, if the adopted value 
— 10 km s^^ were incorrect, the median vr would deviate from zero. The 
root-mean-square scatter of the median vr for subsamples of nearby M stars selected by 
distance and color is 0.5 km s~^, which is an upper limit on the error in the adopted value 
of Vx- This result, which is based on the full three-dimensional velocity distribution, agrees 
well with resu l ts from indirect st atistical deprojection methods using only proper motions 



' Dehnen 


1998; 


Fuchs et al. 


2009) 



For both blue and red disk stars, the extrapolation of the median to Z = 
yields —205 km s~^. Since we corrected stellar velocities for an assumed solar motion of 
—225 km s"^, this implies that the Y component of the solar velocity relative to the bulk 
motion of st ars in the so l ar ne ighborhood is 20 km s~^, in agreement with recent results 



obtained by 



Fuchs et al 



( 120091 ) for the same data set. A similar value was obtained by 
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Dehnen &: Binney ( 1998 ) for theijfl subsample of red stars within 100 pc. 



5. A MODEL FOR THE KINEMATICS OF DISK AND HALO STARS 

Informed by the results from the preceding two sections, we introduce a model that 
aims to describe the global behavior of the observed stellar kinematics. In our model, we 
do not attempt to account for kinematic substructure (e.g., the Monoceros stream), or the 
Galactic bulge region, nor do we incorporate any complex kinematic behavior close to the 
Galactic plane. Nevertheless, we attempt to capture the gross properties of the data in the 
volume probed by SDSS, including the bulk kinematic trends and the kinematic differences 
between high-metallicity disk stars and low-metallicity halo stars. We describe the model 
in § 15. H then test it in § 15.21 using both the proper-motion and radial- velocity samples. 



5.1. The Kinematic Model 

For halo stars, a single velocity ellipsoid (expressed in a Galactocentric spherical 
coordinate system, see § 14. ip is a good description of the gross halo kinematics within the 
10 kpc distance limit of our sample. Our model assumes that the halo has no net rotation 
(see below for a test of this assumption), and that the principal axes of the velocity ellipsoid 
are aligned with a spherical coordinate system. The velocity dispersions measured in §[3] and 
§ in = 135 km s~^, = 85 km s~^, and = 85 km s~^, are expressed in a cylindrical 
coordinate system. The interplay between proper-motion, radial-velocity, and distance 

^They extrapolated the mean azimuthal motion of color-selected samples, which is 
correlated with the radial-velocity dispersion, to zero dispersion and obtained Vy''^'^'^ = 
—5.3 km s~^, used here. 
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measurement errors is complex, so we use Monte Carlo simulations to translate them to 
spherical coordinates; the results of this exercise are = 141 km s^^ and = 75 km s^^ 
{ajjl is unchanged), with uncertainties of ~ 5 km s~^. 

For disk stars within ~ 1 — 2 kpc from the Sun, the velocity-measurement errors 
are sufficiently small, and the samples are sufficiently large, to resolve rich kinematic 
substructure (e.g.. Figure [3]). This behavior is quantified in detail in Kowalski et al. (in 
preparation), but here we simply use the two-component model given by equations f|T9l) 
and ( l20l) to describe the non-Gaussian v^p distribution and velocity shear seen for disk stars. 
Furthermore, we assume that vr and vz have uncorrelated Gaussian distributions, with zero 
mean and the Z-dependent intrinsic dispersion parameters listed in Table [1] As discussed 
in § m there is no compelling evidence for a tilt in the velocity ellipsoid of blue disk stars in 
the vr — Vz plane, so we model the disk velocity ellipsoid in cylindrical coordinates. 



5.2. Global Model Tests 

Our model predicts distributions of the three measured kinematic quantities, frad, f^i, 
and fib, for an arbitrary control volume defined by color, magnitude, and sky coordinates. 
We can test the consistency of this model, m, with our data, d, by computing the residuals 
of each kinematic quantity, 

d — m . . 

where is the measurement error of the data and is the dispersion predicted by 
the model. For all three kinematic quantities, we find x ^ to within 0.05 for all three 
quantities, with dispersions of ~ 1.05 — 1.1. While this result is a necessary condition for 
the model to be acceptable, it is not sufficient. A stronger test of the model was shown in 
Figure [9] for disk stars, but in what follows, we perform additional tests that cover all sky 
regions with available data. 
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5.2.1. Tests With the Radial Velocity Sample 

Figure [16] compares the medians and dispersions for the measured and modeled radial 
velocities of halo stars. The alignment of the velocity ellipsoid with spherical Galactocentric 
radius is clearly seen in the bottom left panel, where the data show a dispersion gradient 
moving away fro m the Galacti c Cent er. The increased dispersion towards / = 180° was 
misinterpreted by llvezic et al.l ( 120061 ) to be a sign of substructure. There are no large 
discrepancies between the measured and predicted behavior - the median value of the 
difference between observed and modeled values is 4.6 km s~^, with a scatter of 19 km s~^ 
(see top right panel). A similar scatter is obtained between two model realizations with the 
same number of stars and measurement errors. The dispersion ratio, o'fad/'^rad; centered 
on 1.13, with a scatter of 0.2 (see bottom right panel). For pairs of model realizations, the 
ratio is always centered on one to within 0.02, with a scatter of 0.2, suggesting that the 
observed velocity dispersion is about 10% larger than predicted by our smooth model. If 
the residuals were due to halo rotation, we would expect a spatial coherence. Similarly, 
for disk stars at 1 < < 2.5 kpc, the median radial velocity residual is 2.8 km s~^, with a 
scatter of 6.6 km s~^ (not shown). 



5.2.2. Tests Based on the Proper Motion Sample 

The large size of the proper-motion sample enables a much higher spatial resolution 
when searching for structure in the model residuals. We have compared the observed 
and modeled proper-motion distributions in narrow bins of distance, across the sky, and 
separately for disk and halo subsamples. As an illustration. Figures [17] and [18] show the 
median longitudinal and latitudinal proper motion observed for disk stars. There is very 
little change in the proper-motion distribution among different distance bins, due to the 
nearly linear vertical rotational- velocity gradient. The residuals for longitudinal proper 
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motion are shown in Figure [191 They provide weak evidence for either substructure or 
a radial gradient that is not modeled, but it is difficult to distinguish between these two 
possibilities with these data alone. 

A comparison of disk and halo subsamples in a distance bin centered on d = A kpc 
is shown in Figure |2D1 The largest data vs. model discrepancy for halo stars, seen in the 
bottom left panel towards Z ~ 0°, is also seen from a different viewing angle in the top 
left panel of Figure [8] (i? ~ 6 kpc and Z ~ 2 kpc) . It is likely that this discrepancy is 
due to contamination of the halo sample by metal-poor disk stars. Figure EH shows the 
residuals for halo stars selected from the 8 — 10 kpc distance bin. The residuals for both 
proper-motion components exhibit similar morphology to the systematic proper-motion 
errors plotted in the two left panels of Figure [TJ In this distance bin, they correspond to 
velocity errors of ~ 30 km s~^; as such, kinematic substructure in the halo will be difficult 
to discern with this sample at > 10 kpc. We note that it is tempting to associate the 
coherent u/ res iduals towards / ~ 300° and I ~ 60° with the Virgo overdensity (see JOS; 



AnetaL 



20091 ). However, the top left panel in Figure [H clearly shows systematic proper 



motion errors of the required amplitude (~ 1.5 mas yr ^) in the same sky region. 



Schlaufman et al. 



(120091 ) use SDSS radial velocity measurements for metal-poor turnoff 
stars to search for pieces of cold debris streams (which they term Elements of Cold Halo 
substructure, or ECHOs). In 137 lines of sight they detect ten ECHOs. The six northern 
ECHOs from their Table 2 (class I peak detections) are shown as white circles in the 
bottom two panels of Figure EH It seems plausible that four of these ECHOs might be 
associated with the Monoceros stream (6 ~ 30°), wh ile one of the two rern a ining detections 



(at / = 162.4° and / = 59.2°) is associated with the iGrillmair fc DionatosI ( l2006l ) stream. 



The ECHO at / = 100.7° and b = 56.8° remains unassociated with any known substructure. 
We conclude that our model reproduces the ffist and second moments of the velocity 
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distributions reasonably well for both disk and halo stars. Except in the region close to the 
Monoceros stream, the non-Gaussian distribution for disk stars is also well-described. 
On average, the model agrees with the data to within ~ 1 mas yr~^ for proper motions, 
and ~ 10 km s~^ for radial velocities. 



5.3. Constraints on i'lsr from Large-Scale Halo Kinematics 

The proper-motion distribution for halo stars towards the Galactic poles depends 
only on the difference between the velocity of the local standard of rest (given that the 
solar peculiar motion is known to ~ 1 km s~^), fLSR^ and for halo stars. At least in 
principle, samples that extend over a large sky area can be used to provide constraints on 
both the halo rotation and f lsr- Figure [22] compares the radial velocity and longitudinal 
proper-motion residuals between two models with (f^^^°,fLSR)= (—20,180) km s~^, and 
(20,220) km s~^ (for both we fixed wlsr ~ = 200 km s~^, to make the models agree 
with the data towards the North Galactic Pole; see Figure [5]). In order to distinguish 
these two models observationally, systematic errors in the radial velocities must be below 
10 km s~^, and systematic errors in the observed proper motions must be below 1 mas yr~^. 
These requirements are comparable to the systematic errors in our data set, so we can only 
state that the data are consistent with wlsr ~ 200 km s~^ (or similarly, no halo rotation) 



to an accuracy q 



(jPerryman et al. 
method. 



20 km s ^ . The proper motion measurements from the Gaia survey 



20011 1 will be sufficiently accurate to exploit the full potential of this 
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5.4. The Kinematic Parallax Relation 

Any constraint on our best-fit model parameters that uses proper motions is sensitive 
to the systematic errors in the distance scale obtained from our photometric-parallax 
relationship. However, since our model is a good fit to the radial velocity data alone (see 
Figure [T6|) . it is possible to estimate distance errors by minimizing the differences between 
the observed and modeled proper-motion distributions. Such a kinematic-parallax relation, 
derived from a combinatio n of radial velocity and prop er-motion data sets, was proposed for 



the solar neighborhood by 



Binney &: Tremaind (119871 ). Unlike in the solar neighborhood. 



our volume is sufficiently large that kinematics vary with position. Nevertheless, it is 
conceptually the same method; a dipole (in our case a more complex angular function) is fit 
to the radial velocity and proper-motion distributions, and the ratio of the best-fit dipole 
magnitudes constrains the distance scale. 

Using only low-metallicity halo stars, we obtain a distance-scale error of ~ 5% - 
our adopted absolute magnitudes should be ~ 0.1 mag brighter. This offset is consistent 
with the expected systematic errors in the calibration of the photometric-parallax relation 
(see 108). In other words, the adopted distance scale properly connects the radial and 
proper-motion distributions. This method provides a much weaker distance-scale constraints 
for disk stars than for halo stars, because the vertical velocity gradient in the disk leads to 
a degeneracy between errors in the distance scale and errors in the adopted disk velocity. 
The radial-velocity measurements do not constrain the velocity scale because the reference 
point depends on distance. 
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SUMMARY AND DISCUSSION 



This is the first analysis based on SDSS data that simultaneously studies the 



kinematics of the halo and disk populations. Past studies of halo stars a. 



performed b y 



Smith et al. 



Sirko et al. 



AUende Prieto et all (120061 ) 



CaroUo et al. 



one were 



(120071 ) and 



(120091 ). while disk samples ran g ing fro r n nearby M stars to dis tant F/G 



stars have been studied by 



Bochanski et al. 



(I2007bh 



West et al. 



(120081), and 



Fuchs et al. 



(120091 ). Throughout this paper, we have quantified the probability distribution function 
P3(/ii, /if,, fradl^ — 9,9,9 — r,l,b), introduced in 108, that describes proper-motion and 
radial- velocity measurements in the 9 vs. 9 — r color-magnitude diagram as a function of 
position on the sky and u — g color. We have developed a simple empirical model with 
disk and halo components that map well to populations detected in the stellar density 



distribution (J08) and the metallicity distribution (108). At distances accessible to the 

lent w ith results from 
(I2OOJ). The extension of 



HIPPARCOS survey (< 100 pc), we ob tain e ncouraging agreement w ith results from 
I I I I 



Nordstrom et al. 



Dehnen fc Binnevi (il998ll . iPehnenl (119981 ). and 
kinematic mapping to distances up to ~ 10 kpc with millions of stars represents a significant 
observational advance, and delivers powerful new constraints on the dynamical structure of 
the Galaxy. In less than two decades, the observational material for such in situ mapping 



has p rogressed from first pioneerin g studies based on on ly a few hundred objects (IMajewski 



19921 ). to over a thousand objects (IChiba fc Beers! l2000l ). to the massive data set discussed 
here. 



6.1. Disk Kinematics 



The disk kinematics of the Milky Way are dominated by rotation with a s mooth vertica l 
gradient. Our analysis extends the e arly measurements of this gradient (e.g., lMurraylll986 : 
Majewskilll992l : IChiba fc Beersll2000l ) to vastly larger sky area and to a much larger distance 
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range. The mean rotational velocity and the three velocity dispersions for disk stars can be 
modeled as simple functions of the form a + h\Z\^ (see equations f[T^ and f l2U]) . and Tabled]). 
The rotational velocity distribution for the disk component is non-Gaussian, and can be 
formally modeled as a sum of two Gaussian components with fixed normalization ratio for 
0.1 < |Z|/kpc < 4. The fact that the normalization ratio of these two components does 
not vary with Z is at odds with the standard disk decomposition into thin- and thick-disk 
components ( see also sections 3.4 . 4 and 4.2.1 in 108). B ased on N-body simulations 



performed by 



Roskar et al 



mm 



Loebman et al. 



(|2008[ ) argued that the absence of a 



velocity- metallicity correlation at the thin/thick disk boundary, pointed out by 108, may be 



due 
also 



;o a combination of a stro ng vertical age gradient and the radial migration of stars (see 



Schonrich fc Binney 



20091 ). A more detailed study will be presented in Loebman et al. 
in preparation). A significant vertical a ge gradient for d i sk stars is also supported by an 



analysis of active M dwarfs presented in iBochanski et al.l (l2007b( ). Such an age gradient, 
together with the measured velocity dispersion-age correlations for local disk stars (e.g.. 



Nordstrom et al. 



2004; 



Rocha-Pinto et al. 



2004 



West et al. 



20081 ). may be responsible for 



the measured increase of velocity dispersions with distance from the Galactic plane. 

Close to the plane, the proper-motion data imply a complex multi-modal velocity 
distribution that is inconsistent with a description based on a simple Schwarzschild ellipsoid. 
It is reassuri ng that we ob tained a velocity-distribution morphology very similar to that 



obtained by 
directly by 



Dehnen ( 



19981) using s tatistical deprojection of the HIPPARCOS data, and 



Nordstrom et al. 



( 12004J ) using a subsample of HIPPARCOS stars with full 
three-dimensional velocity information. In addition, our results for the first and second 
moments of the velo city distribution fo r nearby M stars agree with analogous results 



obtained recently by 



Fuchs et al 



(12009[ ). The orientation of the velocity ellipsoid is strongly 
affected by multi-modal structure, so one should take care in its interpretation. We discuss 
these issues in more detail in Kowalski et al. (in preparation). 



-44- 



6.2. Halo Kinematics 



Our results fo r the velocity distribution of halo stars are in excellent agreement with 



Smith et al. 



(|2009[ ). They used proper-motion measurements based only on SDSS data, 
and thus have significantly different, and most likely much smaller, systematic errors 
than t he SDSS-PO S S pro per-motion measurements analyzed here. The much larger size 



of the 



Munn et al. 



( I2OO4J ) catalog analyzed here allows us to rule out the possibility 
that the Smith et al. result was biased by local substructure. The close agreement 
of our results for the orientation and size of the halo velocity ellipsoid (we obtained 
= 141 km s~^, = 75 km s^^, and (7¥ = 85 km s~^, while thei r values are 142 krn s~^ 



Carollo et al. 



2OO9I). 



77 km s^^, and 81 km s~^, respectively) are encouraging (see also 
Their estimated errors of 2 km s^^ apparently do not include systematic effects (such as 
errors in photometric parallax; both studies used the same calibration from 108) - based 
on our Monte Carlo simulations, we believe that the true errors cannot be smaller than 
~ 5 km s^^. Additional independent evide nce for the tilt of the velocity ellipsoid comes 



from the RAVE survey (ISiebert et al 



20081 ) . These measurements of the velocity ellipsoid 



for halo stars rep resent a strong constra i nt for the shape of gravi t ationa l pot ential, as 



discus sed by, e.g.. 



fcoooh . 



Amendt fc Cuddefordl (|l99lh . 



Kent fc de Zeeuwl (jl99l[ l and 



Smith et al. 



We note that 



Majewski 



(119921) measured a retrograde halo rotation using stars observed 
towards the North Galactic Pole (in our nomenclature, he obtained a mean rotational 
velocity w$ = 50 ± 16 km s~^). A star- by-star comparison of his data and the data analyzed 
here showed that photometric and proper-motion measurements agree within the stated 
errors. The main reason for different conclusions about halo rotation are the different 
distance scales: his distances ar e on average 30% la rger than our distances, resulting in 



larger tangential velocities. The 



Carollo et al. 



(120071 ) claim for a large outer-halo retrograde 
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rotation, further refined by ICaroUo et al.l ( 120091 ) (in our nomenclature, a mean rotational 
velocity f$ = 80 ± 13 km s^^), remains intact, as their distance scale is similar to ours. 

The kinematic measurements for halo stars presented here should not be extrapolated 



beyond the sample distance limit of 10 kpc. For example, using 24 



stars, globular clusters, and satellite galaxies. 



Battaglia et al. 



halo objects, including 



( I2OO5I ) detected a continuous 



decline of the radial velocity dispersion beyond a Galactocentric radius of ~ 30 kpc, from 
about 120 km s~^ to 50 km s~^ at ~ 120 kpc. In addition, the distance limit of our 
sample, together with the decreasing sensitivity of the photometric metallicity indicator for 
[Fe/H] < —2. 5 , prevent u s from robustly testing the possible halo dichotomy discussed by 



CaroUo et al. 



(2007 



20091 ). 



6.3. Kinematic Substructure 



The model developed h ere can be used to search for kinematic substructure with a low 



contrast level. For example, iSchlaufman et al. 



(120091 ) had to generate a background model 



when searching for cold streams 
neighborhood was performed by 



using radial ve . 



Klement et al. 



ocity data; a similar study for the solar 



(I2OO9I ). A global model-based description is 



especially important when using large numbers of proper-motion measurements to search 
for substructure. While radial-velocity data are superior at large distances, searches based 
on proper motions should be better within a few kpc, due to the high completeness and 
much larger sample size. A user-friendly interface to our model code Galfasu, which allows 
generation of mock catalogs in an arbitrary direction (or across the entire sky) and to an 
arbitrary depth, will be described elsewhere (Juric et al., in preparation). 



The Monoceros stream is clearly detected as a major outlier from the smooth model 



^Please see www.mwscience.net/galfast 
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presented here. We have not found any other large kinematic substructure within 10 kpc 
that deviates at a detectable level. Other deviations from the smooth model predictions 
are likely due to systematic errors in the proper motion and radial velocity measurements. 
However, since the main goal of this paper was to quantify the overall kinematic behavior, 
the emphasis of our analysis was on the first and second moments of the kinematic 
qu antities. It is hkely that mor e sophisticated s tatist ical methods, such as those discussed 



by 



Schlaufman et al. 



(120091 ) and 



Klement et al 



(|2009[ ). will be more efficient in se arching for 



substructure. As an example, we have verified that a moving group discovered by 



Majewski 



( 119921 ) is reproducible with our data. 



6.4. Future Surveys 



The results presented here will be greatly extended by several upco ming large-scale. 



deep optical 



survey s, including the Dark Energy Survey (IFlaughei 



( iKaiser et al. 



20021 ). and the Large Synoptic Survey Telescope (llvezic et al. 



2008) 



'an-STARRS 



2008bh . These 



surveys will extend the faint limit of this sam ple and that of the upcoming Gaia mission 



( iPerryman et al 



2001 



Wilkinson et al. 



20051 ) by 4 — 6 mag. For example, LSST will obtain 



proper motion measurements of comparable accuracy to those of Gaia at their faint limit, 
and smoothly extend the error vs. magnitude curve deeper by 5 mag (for details see Eyer 
et al., in preparation). With its u-band data, LSST will enable studies of metallicity 
and kinematics using the same stars out to a distance of ~ 100 kpc (~ 200 million F/G 
main sequence stars brighter than g = 23.5, for a discussion see 108). By comparison, the 



best measurements of the outer- 



several hundred ( iBattaglia et al. 



lalo radial velocity dispers ion to date are based on from 



120051 ) to several thousand ( IXue et al.ll2008l ) objects. These 
upcoming studies are thus certain to provide valuable new information about the formation 
and evolution of our Galaxy. 
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A. The Revised SDSS Metallicity Scale 



Analysis of the metallicity and kinematics of halo and disk stars by 108 utilized 
photometric metallicity estimates for F/G stars with 0.2 < g — r < 0.6. Their mapping 
function from the g — r vs. u — g color-color diagram to metallicity was calibrated using 
stars with spectroscopic metallicities distributed in SDSS Data Release 6. At that time. 



high-metallicity stars required for the calibration of methods implementec 



Beers et al 



in the automated 



20061) were not 



spectroscopic pipeline (SEGUE Stellar Parameters Pipeline, 
available. Between Data Releases 6 and 7, the required data were collected and the new 
calibration resulted in the re vised spectroscopic metalli city values distributed with Data 



Release 7 ( iLee et al. 



2008a 



b 



AUende Prieto et al. 



20081 1. 



Here, we recalibrate the photometric-metallicity estimator using updated spectroscopic 
metallicities from Data Release 7. In addition, we re-derive the parts of the 108 analysis 
that are most affected by this change in the metallicity scale. 



A.l. The Updated Photometric-Metallicity Estimator 

As shown in Figure lAlt the largest difference between the SDSS spectroscopic 
metallicity values distributed with Data Releases 6 and 7 is, as expected, at the 
high-metallicity end. In particular, the abrupt cutoff in the metallicity distribution at 
[Fe/H] ~ —0.5 (see figure 9 in 108) is no longer present, and the distribution extends to 
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values as high as [Fe/H] ~ —0.2 (the distances for the stars shown range from ~ 1 kpc to 
~ 7 kpc). 

We proceed to re-derive the photometric-metalhcity cahbration using the same selection 
criteria and the same methodology as in 108. The new data set admits a slightly simpler 
function - the double definition of the x axis is no longer required, and the new expression 
is 

[Fe/H]ph = A + Bx + Cy + Dxy + Ex^ + Fy"^ + Gx^y + Hxy"^ + Ix^ + J?/^ (Al) 

with X = {u — g) and y = [g — r). The best-fit coefficients are {A — J) = (—13.13, 14.09, 
28.04, -5.51, -5.90, -58.68, 9.14, -20.61, 0.0, 58.20). Note that the coefficient for x^, I, is 
to zero. We removed this term because it was producing too much curvature at the right 
end (red u — g) of the best-fit map. 

We estimate that an upper limit on the intrinsic metallicity scatter for fixed noiseless 
u — g and g — r colors (presumably due to limited sensitivity of broadband colors to 
metallicity variations) is about 0.1 dex. This value is estimated from the scatter in the 
difference between spectroscopic and photometric metallicities, discussed below. Unlike 108, 
who simply adopted the median metallicity value given by the above expression for each 
star, we draw photometric metallicity estimates from a Gaussian distribution centered on 
the best-fit median value, and with a width of 0.1 dex. The main benefit to this is that we 
avoid hard edges in the photometric metallicity distribution for stars close to the edges of 
the calibration region in the g — r vs. u — g diagram. 

The performance of the new map is qualitatively similar to that of the old map. The 
mean and rms values for the difference between spectroscopic and photometric metallicities 
as functions of the g — r and u — g colors are shown in the top two panels of Figure IA2[ 
Typical systematic errors in the map (i.e. the median difference per pixel) are ~ 0.1 dex or 
smaller, and the scatter varies from ~ 0.2 dex at the high-metallicity end to ~ 0.3 dex at 



- 50 - 



the low-metallicity end (note that this scatter includes a contribution from errors in both 
spectroscopic and photometric metalhcities). 

The above photometric metalhcity estimator is apphcable for stars with 0.2 < g—r < 0.6 
and —0.25 + 0.5(m — g) < g — r < 0.05 + 0.5('U — g); that is, for main-sequence F and G 
stars in the cahbration region of the g — r vs. u — g color-color diagram (top two panels 
of Figure E2|) . For stars with spectroscopic metallicity [Fe/H] > —2.2, the distribution of 
the difference between spectroscopic and photometric metalhcities is well- described by a 
Gaussian with a width of 0.26 dex (see the bottom right panel in Figure lA2|) . 

It should be noted that the performance of the photometric-metallicity estimator 
deteriorates at the low-metallicity end because the u — g color becomes insensitive to 
decreases in metallicity. As shown in the bottom left panel of Figure \A2\ the photometric 
metallicity saturates at [Fe/H] ~ —2 for smaller values of spectroscopic metallicity. Even 
at [Fe/H] = —2, the true metallicity is overestimated by 0.2 — 0.3 dex, and by [Fe/H] = —3 
this bias is as large as 1 dex (the photometric-metallicity values never become significantly 
lower than [Fe/H] = —2). This shortcoming could be partially alleviated by employing more 
accurate w-band photometry (say, with errors of 0.01 mag instead of 0.03 mag, as used here), 
but probably not for metalhcities lower than [Fe/H] = —2.5. Fortunately, the low-metallicity 
inner-halo stars within SDSS reach have a median metallicity of [Fe/H] ~ —1.5 (108), 
so for the majority of stars the photometric metalhcities are robust. Another important 
note is that, despite the improvement at the high- metallicity end, the calibration range 
only extends to [Fe/H] ~ —0.2. Our calibration sample did not include young, metal-rich 
main-sequence stars {u — g > 1.1 and g — r < 0.3). For this reason, our polynomial model 
underestimates true metalhcities by 0.2 — 0.3 dex at [Fe/H] = 0, and probably more for 
[Fe/H] > 0. Any result relying on higher metalhcities should be interpreted with caution 
(especially at low Galactic latitudes where the uncertain ISM extinction may strongly affect 
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the estimated metallicities). We plan to extend our calibration f urther into the met al-rich 



domain by employing data from the SDSS-III SEGUE-2 survey (lYanny et al. 



20091), and 



further refinements of the SEGUE Stellar Parameters Pipeline, which are now underway. 



A. 2. Tomography II Reloaded 

108 pointed out several aspects of their analysis that may have been affected by the 
metallicity "compression" at the high-metallicity end in DR6. We repeated their full 
analysis and report here on those aspects where differences warrant discussion. 

The "hard" upper limit on photometric metallicity estimates at the high-metallicity 
end ([Fe/H] ~ —0.5) with the DR6 calibration is best seen in the bottom left panel in 
figure 9 from 108. We reproduce that map of the conditional metallicity distribution in the 
top left panel of Figure IA3I As expected, the metallicity distribution of disk stars within 
2 kpc of the Galactic plane now extends to [Fe/H] ~ 0. 

In the new calibration, the parameters describing the variation of the median metallicity 
for disk stars as a function of the distance from the Galactic plane, 

fioiZ) = /ioo + exp{-\Z\/H^) dex, (A2) 

are also changed. The updated values are = 0.5 kpc, /ioo = —0.82 and = 0.55 (the 
old values were = 1.0 kpc, /ioo = —0.78 and = 0.35). The best-fit values of /ioo 
and A^ are accurate to ~ 0.05 dex. Values of in the range 350 — 700 pc are consistent 
with the data - the decrease in Hf, is required by the local constraint, Hd^Z) = —0.2 



( iNordstrom et al. 



2004 



Allende Prieto et al. 



20041 ). 



An interesting result from 108 was the detection of disk stars at a distance from the 
Galactic plane as large as ~ 6 kpc (see their Figure 10). A peak at [Fe/H] = —0.5 in 
the metallicity distribution of stars at those distances was another manifestation of the 
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metallicity "compression." As demonstrated in the top right panel of Figure IA3t this 
peak is not present when using the revised cahbration. However, there is still statistical 
evidence that disk stars exist at these large distances from the plane: about 5% of stars in 
the 5 < Z < 7 kpc bin are consistent with disk sta rs, in agreem e nt wi th ex trapolation of 



the ex ponential profile derived from stellar counts. 



Lee fc BeersI (120091 ) and 



CaroUo et al. 



( I2OO9I ) have also commented on signatures of thin-disk-like chemistry and kinematics for a 



small fraction of stars several kpc above the Galactic plane. 

Perhaps the most intriguing result of 108 was the non-detection of a correlation between 
rotational velocity and metallicity for disk stars at Z ~ 1 kpc. At such distances from 
the Galactic plane, the counts of thin- and thick-disk stars are expected to be similar. 
Since traditionally the thick disk component is associated with a larger velocity lag and 
lower metallicities, a fairly strong and detectable correlation was expected (see 108 for 
details). The two bottom panels in Figure 1X3] demonstrate that such a correlation is still 
undetected, although the photometric metallicity range now extends to higher values (up 
to [Fe/H] ~ -0.2). 

The higher metallicity values obtained with the re-calibrated relation have a 
quantitative effect on the best-fit metallicity distributions shown in Figure 7 of 108. Using 
the same methodology, we reproduce the metallicity distributions with the new calibration 
in Figure IA4[ 108 modeled the non-Gaussian disk metallicity distribution using a sum of 
two Gaussians with a fixed amplitude ratio (1.7:1), fixed difference of the mean values 
(0.14 dex), and fixed widths (0.21 dex and 0.11 dex), which "slides" as a function of Z, 
according to equation (1A2|) . We find that the only significant change is an increase to 
the width of the second Gaussian to 0.21 dex, which accounts for the extension of the 
metallicity distribution to higher values. Only minor changes are required for the best-fit 
halo metallicity distribution (see Table 3 in 108): the median halo metallicity is now —1.46 
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in the first three bins, and —1.56 in the most distant Z bin, and its width changed from 
0.32 dex to 0.36 dex in the last bin. We note somewhat less scatter of the data points 
around the best-fit functions with the re-calibrated data set. To summarize, the revised 
best-fit parameters that describe halo and disk metallicity distributions are: 

• The halo metallicity distribution is spatially invariant and well-described by a 
Gaussian distribution centered on [Fe/H] = —1.46, with a width an = 0.30 dex 
(not including measurement errors). For \Z\ < 10 kpc, an upper limit on the halo 
radial-metallicity gradient is 0.005 dex kpc~^. 

• The disk metallicity distribution varies with Z such that its shape remains fixed, 
while its median, fi£,, varies as given by equation flA2|) (with best-fit parameter values 
Hfj^ = 0.5 kpc, fioo = —0.82, and = 0.55). The shape of the disk metallicity 
distribution can be modeled as 

poix = [Fe/H]|Z) = 0.63G[x|/i = a{Z),a = 0.2] + 0.37 = a{Z) + O.U,a = 0.2], 

(A3) 

where the position a and the median fj,D are related via a{Z) = iid{Z) — 0.067 (unless 
measurement errors are very large). 

We point out that the asymmetry of the metallicity distribution for disk stars 
is now less pronounced (as implied by the same widths of the two best-fit Gaussian 
components). Nevertheless, due to our large sample size, the non-Gaussianity is detected 
at high significance. A remaining uncertainty is the error distribution of the photometric 
metallicities, which itself could account for such a deviation from Gaussianity. However, 
to the extent possible using a highly incomplete spectroscopic sample (c.f. the bottom 
right panel in Figure R2] and discussion in 108), we are unable to quantitatively explain the 
observed deviation from Gaussianity as an artifact of the photometric-metallicity method. 
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B. SQL Query Example 

The following SQL query was used to select and download data for all SDSS stars with 
spectroscopic and proper- motion measurements (see http://casjobs.sdss.org/CasJobs). 

SELECT 

round(p.ra,6) as ra, round (p. dec, 6) as dec, 

p. run, p. rerun, round(p.extinction_r,3) as rExt, 

round(p.psfMag_u,2) as upsf , comments are preceded by 

round(p.psfMag_g,2) as gpsf , rounding up 

round(p.psfMag_r,2) as rpsf, 
round(p.psfMag_i,2) as ipsf, 
round(p.psfMag_z,2) as zpsf, 
round(p.psfMagErr_u,2) as uErr, 
round(p.psfMagErr_g,2) as gErr, 
roimd(p.psfMagErr_r,2) as rErr, 
round(p.psfMagErr_i,2) as iErr, 
round(p.psfMagErr_z,2) as zErr, 

round(s .pmL,2) as pmL, round(s .pmB,2) as pmB, pmRaErr, 
t.specObjID, t. plate, t.mjd, t.fiberid, t.feha, t.fehaerr, 
t.logga, t.loggaerr, t.elodierv, t . elodierverr 
INTO mydb.pmSpec 

FROM star p, propermotions s, sppParams t, specobjall q 
WHERE 

p.objID = s.objID and s. match = 1 must have proper motion 

and t.specobjid = q.specobjid and q.bestobjid = p.objid 
and s.sigra < 350 and s.sigdec < 350 quality cut on pm 
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and (p. flags & '4295229440') = see text for flag list 

and p.psfMag_r > 14.5 avoid saturation 

and p.psfMag_r < 20 practical faint limit for pm 

the end of query 

The following SQL query was used to select and download data for all spectroscopically 
confirmed quasars with proper-motion measurements and redshifts in the range 0.5 — 2.5. 

SELECT 

round (p. ra, 6) as ra, round (p. dec, 6) as dec, 
p. run, p. rerun, round(p.extinction_r,3) as rExt, 
round(p.psfMag_u,2) as upsf, 
round(p.psfMag_g,2) as gpsf, 
round(p.psfMag_r,2) as rpsf, 
round(p.psfMag_i,2) as ipsf, 
round(p.psfMagErr_u,2) as uErr, 
round(p.psfMagErr_g,2) as gErr, 
round(p.psfMagErr_r,2) as rErr, 
round(p.psfMagErr_i,2) as iErr, 

round(s.pmL,2) as pmL, round(s .pinB,2) as pmB, pmRaErr, 
q.specObjID, q. plate, q.mjd, q.fiberlD, 
q.z, q.zErr, q.zConf, q.zWarning, q.specClass 
INTO mydb.pmQSO 

FROM star p, propermotions s, specobjall q 
WHERE 

p.objID = s.objID and s. match = 1 
and q.bestobjid = p.objid 
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and s . sigra < 350 and s . sigdec < 350 — per \citet{Munn04} 
and (p. flags & '4295229440') = 
and p.psfMag_r > 14.5 
and p.psfMag_r < 20 

and q.z > 0.5 redshift limits 

and q.z < 2.5 
the end of query 
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Table 1. Best-Fit Parameters" for the Disk Velocity 
Distribution^. 



Quantity 


a 


h 


c 




-194 


19.2 


1.25 




12 


1.8 


2.0 




34 


1.2 


2.0 




30 


3.0 


2.0 




40 


5.0 


1.5 




25 


4.0 


1.5 



The uncertainties are typically ~ 10 km s ^ for a, ~20% for b and 0.1-0.2 for c. 

^AU hsted quantities are modeled as a + h\Z\'^, with Z in kpc, 
and velocities in km s~^. 

''The distribution is non-Gaussian, and can be formally de- 
scribed by a sum of two Gaussians with a fixed normalization ratio 
/fc:l, with fk — 3.0. The mean value for the second Gaussian has 
a fixed offset from the first Gaussian, v^'^ — {H^^ — A'i^), with 
= 34 km/s. Extrapolation beyond Z > b kpc is not reliable. 
The velocity dispersion for the second Gaussian is given by cr^. If 
this non-Gaussianity is ignored, the v^i^ dispersion is given by erf*. 
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Table 2. A Model Prediction for Velocity Ellipsoid Tilf^ 



R (kpc) 


Z (kpc) 


Orz (deg) 


arctan(Z//2) (deg) 


8.0 


0.7 


3 


5.0 


8.0 


2.0 


11 


14.0 


8.0 


3.5 


26 


23.6 


7.0 


4.0 


32 


29.7 


9.0 


4.0 


25 


24.0 


11.0 


4.0 


21 


20.0 



^The first two columns determine the position in the Galaxy. 
The third column lists the predicted orientation of the velocity 
ellipsoid computed as described in § I4.1.1[ The fourth column lists 
the orientation of a velocity ellipsoid pointed towards the Galactic 
center. 




Fig. 1. — Behavior of proper-motion measurements for 60, 000 spectroscopically confirmed 
SDSS quasars with 6 > 0°. The color-coded maps (see the legend on top, units are mas yr^^) 
show the distribution of the median (left) and rms (right) for the longitudinal (top) and 
latitudinal (bottom) proper motion components in a Lambert projection of the northern 
Galactic cap. The median number of quasars per pixel is ~ 250. For both components, the 
scatter across the sky is 0.60 mas yr~^. The median proper motion for the full quasar sample 
is 0.15 mas yr^^ in the longitudinal direction, and —0.20 mas yr^^ in the latitudinal direction. 
The thick line in the top left panel shows the selection boundary for the "meridional plane" 
sample. 




Fig. 2. — Change of the w$ vs. vr velocity distribution with metalhcity, at an approximately 
constant R and Z. Velocities are determined from proper-motion measurements. The top 
left panel shows the f$ vs. vr diagram for ~ 6, 000 blue (0.2 < g — r < 0.4) stars at 
Z = A — 5 kpc and detected towards the North Galactic Pole {b > 80°). The distribution is 
shown using linearly-spaced contours, and with color-coded maps showing smoothed counts 
in pixels (low-to-high from blue-to-red). The other five panels are analogous, and show 
subsamples selected by metallicity, with the [Fe/H] range listed above each panel (also 
listed are the median r-band magnitude and subsample size). The measurement errors are 
typically 70 km s~^ (per star). Note the strong variation of median f$ with metallicity. 




Fig. 3. — Similar to Figure [21 except that the vs. vr velocity distribution is plotted for 
a range of Z. The top row shows the vq, vs. vr diagrams for ~ 60, 000 red {g — r > 0.6) 
stars at Z = 100 — 700 pc and observed towards the North Galactic Pole. Each panel 
corresponds to a narrow range in Z, given above each panel. The measurement errors 
vary from ~ 3 km s~^ in the closest bin to ~ 12 km s~^ in the most distant bin. Note 
the complex multi-modal substructure in the top left panel. The bottom three panels are 
analogous, and show the f$ vs. vr diagrams for ~ 7, 000 blue (0.2 < g — r < 0.4) stars with 
high metallicity {[Fe/H] > —0.9). The measurement errors vary from ~ 20 km s~^in the 
closest bin to ~ 35 km s~^in the most distant bin. Note that the median f $ approaches zero 
as Z increases. 
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Fig. 4. — Similar to the top left panel in Figure [H except that stars are selected from a 
distance bin that correspo nds to HIPPA RCOS sample (Z = 50 — 100 pc). The positions 



of Eggen's moving groups (lEggen 



19961 ) are marked by circles, according to the legend in 
the bottom right corner. The horizontal line at Vs = —225 km s~^ corresponds to vanishing 



heliocentric motion in the rotational direction. 
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Fig. 5. — Dependence of the rotational velocity, f<j,, on distance from the Galactic plane 
for 14,000 high-metallicity {[Fe/H] > -0.9; top left panel) and 23,000 low-metalhcity 
{[Fe/H] < —1.1, top right) stars with b > 80°. In the top two panels, individual stars are 
plotted as small dots, and the medians in bins of Z are plotted as large symbols. The 2a 
envelope around the medians is shown by dashed lines. The bottom two panels compare the 
medians (left) and dispersions (right) for the two subsamples shown in the top panels and 
the dashed lines in the bottom two panels show predictions of a kinematic model described 
in text. The dotted lines in the bottom right panel show model dispersions without a 
correction for measurement errors (see Table [1]) . 
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Fig. 6. — Symbols with error bars are the measured rotational- velocity distribution, 
for stars with 0.2 < ^ - r < 0.4, h > 80°, and Z = 0.8 - 1.2 kpc (top left, ~ 1,500 
stars), 1.5 — 2.0 kpc (top right, ~ 4, 100 stars), 3.0 — 4.0 kpc (bottom left, ~ 6,400 stars) and 
5.0 — 7.0 kpc (bottom right, ~ 12, 500 stars). The red and green curves show the contribution 
of a two-component disk model (see Equat ions fT9l and 1201) . the blue curves show the Gaussian 
halo contribution, and the magenta curves are their sum. Note the difference in the scale of 
the y axis between the top two panels. 
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Fig. 7. — Analogous to Figure O but for the radial- velocity component, vr. 
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Fig. 8. — Dependence of velocity, measured using proper motions, on cylindrical Galactocen- 
tric coordinates for 172, 000 metal-poor halo stars {[Fe/H] < —1.1; top panels) and 205, 000 
metal-rich disk-like stars {[Fe/H] > —0.9; bottom panels). Stars are selected from three 
regions: b > 80° (the North Galactic pole), 170° < / < 190° (anticenter), and 350° < / < 10° 
(Galactic center). The left column plots rotational velocity, while the right column plots 
vb = sin{b)vii + cos{b)vz- To aid visualization of these boundaries, see the thick line in 
the top left panel in Figure [1] The median values of velocity in each bin are color-coded 
according to the legend shown in each panel. The measurements are reliable for distances 
up to about 7 kpc, but regions beyond this limit are shown for halo stars for completeness. 
The fraction of disk stars is negligible at such distances; their velocity distribution is shown 
for Z < 6 kpc. The region with negative velocity on the right side of top left panel is due 
to contamination of the halo sample by stars from the Monoceros stream. The thin region 
with negative velocity on the left side of top right panel is a data artifact (see text). 




Fig. 9. — Left panel is analogous to the bottom left panel in Figure [H but for the model 
described in the text. The right panel shows the median difference between the data and 
model. Large discrepancies at i? > 12 kpc are due to the Monoceros stream (at R = 18 kpc 
and Z = 4 kpc; disk stars rotate with a median ~ —100 km s~^, while for the Monoceros 
stream, v^f, ~ —200 km s~^). 
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Fig. 10. — Distribution of stars with 0.2 < g — r < 0.4 and Z = 4 — 6 kpc in the rotational 
velocity vs. metallicity plane, for four ranges of Galactocentric cylindrical radius, R (top left: 
3 — 4 kpc; top right: 7 — 9 kpc; bottom left: 12 — 13 kpc; bottom right: 17— 19 kpc). In each 
panel, the color-coded map shows the logarithm of counts in each pixel, scaled by the total 
number of stars. The horizontal lines at f $ = km s~^ and f$ = —220 km s~^ are added to 
guide the eye. High-metallicity {[Fe/H] ~ —1) stars with fast rotation (fij, ~ —220 km s^^) 
visible in the bottom right panel belong to the Monoceros stream, and are responsible for 
the features seen at i? > 15 kpc in the two left panels in Figure [H 
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Fig. 11. — Similar to Figure |5l but for the vertical- velocity component, vz, and using a 
sample of stars with SDSS radial-velocity measurements, 0.2 < g — r < 0.4 and b > 20° 
(12, 000 stars in the high-metallicity subsample, and 38, 000 stars in the low-metallicity 
subsample). An analogous figure for extended samples of 53, 000 disk stars and 47, 000 halo 
stars with 0.2 < g ~ r < 0.6 has a similar appearance. The behavior of the rotational- and 
radial- velocity components in this sample is consistent with that shown in Figures [5] and [71 




Fig. 12. — Three two-dimensional projections of the velocity distribution for two subsamples 
of candidate halo stars selected using spectroscopic metalhcity (—3 < [Fe/H] < —1.1) and 
with 6 < -R/kpc < 11. The top row corresponds to 2, 700 stars with distances, 3 < Z/kpc < 
4, and the bottom row to 1,300 stars with —4 < Z/kpc < —3. The distributions are 
shown using linearly-spaced contours, and with a color-coded map showing smoothed counts 
in pixels (low-to-high from blue-to- red) . The measurement errors are typically 60 km s~^. 
Note the strong evidence for a velocity-ellipsoid tilt in the top and bottom left panels (see 
also Fig. [13]). The two dashed lines in these panels show the median direction towards the 
Galactic center. 




Fig. 13. — Illustration of the velocity-ellipsoid tilt- angle variation. Analogous to Fig. [ 
except that only the Vz vs. Vr projection is shown for a constant Z, for three ranges of 
as marked on the top of each panel. 




Fig. 14. — Analogous to Figure [121 except that the velocity distribution is shown for two 
subsamples of candidate disk stars selected using spectroscopic metallicity {[Fe/H] > —0.9). 
The top row corresponds to 2, 200 stars with distances from the Galactic plane 1.5 < Z/kpc < 
2.5, and the bottom row to 1, 500 stars with —2.5 < Z/kpc < —1.5. The measurement errors 
are typically 35 km s~^. Note the absence of velocity-ellipsoid tilt in the top and bottom left 
panels. 




Fig. 15. — Analogous to Figure [HI except that the velocity distribution is shown for two 
subsamples of red stars {g — r > 0.6): the top row corresponds to 3, 000 stars with distances 
from the Galactic plane 0.6 < Z/kpc < 0.8, and the bottom row to 4, 600 stars with —0.8 < 
Z/kpc < —0.6. The measurement errors are typically ~ 15 km s~^. There is no strong 
evidence in these panels for a tilt in the velocity ellipsoid. 
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Fig. 16. — Comparison of medians and dispersions for the measured and modeled radial 
velocities of 20,000 blue (0.2 < g — r < 0.4) halo stars (spectroscopic [Fe/H] < —1.1) at 
distances, D = 2 — 7 kpc, and b > 20°. The top left panel shows the median measured 
radial velocity in each pixel, color-coded according to the legend shown at the top (units are 
km s~^). The top right panel shows the difference between this map and an analogous map 
based on model-generated values of radial velocity, using the same scale as in the top left 
panel. The bottom left panel shows the dispersion of measured radial velocities, color-coded 
according to the legend above it. The bottom right panel shows the ratio of this map and 
an analogous map based on model-generated values of radial velocity, color-coded according 
to the legend above it. When the sample is divided into 1 kpc distance shells, the behavior 
is similar. 




Fig. 17. — Distribution of the median longitudinal proper motion in a Lambert projection of 
the North Galactic cap for high-metallicity {[Fe/H] > —0.9), blue (0.2 < g — r < 0.4) stars, 
in several distance bins (top left: 58,000 stars with D = 0.8 — 1.2 kpc; top right: 119,000 
stars with D = 1.8 - 2.2 kpc; bottom left: 72, 000 stars with D = 2.8 - 3.2 kpc; bottom 
right: 43, 000 stars with D = 3.8 — 4.2 kpc). All maps are color-coded using the same scale, 
shown in the middle (units are mas yr~^). Note that the magnitude of the proper motion 
does not change appreciably as the distance varies from ~ 2 kpc to ~ 4 kpc; this is due to 
a vertical gradient in the rotational velocity for disk stars (see Figure |3]). 
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Fig. 18.— 



Analogous to Figure HTJ except that the latitudinal proper motion is shown. 




Fig. 19. — Similar to Figure [T71 except that the median difference between the observed 
value of longitudinal proper motion and a value predicted by the model described in the text 
is shown. All maps are color-coded using the same scale, shown in the middle. Note that 
the displayed scale is stretched by a factor of two compared to the scale from Figure [T71 in 
order to emphasize discrepancies. 




Fig. 20. — Similar to Figures [T71 and [T9| except that the behavior of high-metalhcity (left) 
and low-metallicity (right) stars is compared in a single distance bin (3.5 — 4.5 kpc). The 
top two panels show the median longitudinal proper motion, and the two bottom panels 
show the median difference between the observed and model-predicted values. An analogous 
figure for the latitudinal proper motion has similar characteristics. 




Fig. 21. — Similar to Figure |20| except that the behavior of low-met alhcity stars from the 
8 — 10 kpc distance bin is analyzed. The top two panels show the median longitudinal (left) 
and latitudinal (right) proper motions, and the two bottom panels show the median difference 
between the observed and model-predicted values. The maps are color-coded according to the 
legends in the middle (mas yr~^; note that the bottom scale has a harder stretch to emphasize 
structure in the residual maps). The two bottom panels display very similar morphology to 
systematic proper-motion errors shown in the two left panels in Figure [H In the bottom 
panels, the white symbols show the positions of the six northern cold substructures identified 
by Schlaufman et al. (2009). 




Fig. 22. — Comparison of radial- velocity residuals (left panel; analogous to the top right 
panel in Fig. [16]) and longitudinal proper-motion residuals (right panel; analogous to the 
bottom left panel in Fig. [201 except for the larger distance range) for two halo models with 
(^^J"'°,t^LSR)= (-20,180) km s-\ and (20,220) km s-\ Note that we set vlsr - vf^° = 
200 km s~^. The residuals are color-coded according to the legend above each panel (units 
are km s~^ for the left panel and mas yr~^ for the right panel). In order to distinguish these 
models, systematic errors in radial velocity must be below 10 km s^^, and systematic errors 
in proper motion must be below 1 mas yr~^. 
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Fig. Al. — Summary of the differences between SDSS spectroscopic-metallicity values dis- 
tributed with Data Releases 6 and 7. The left panel shows the median difference between 
the DR7 and DR6 values for 0.02 x 0.02 mag^ bins in the g — r vs. u — g color-color di- 
agram, color-coded according to the legend shown in the panel. The largest differences of 
0.2 — 0.3 dex are seen in the top right corner, which corresponds to high metallicities. The 
right panel shows the difference in metallicities as a function of the new DR7 values. Indi- 
vidual stars are shown as small dots, and the median values of the difference are shown as 
large circles. The two dashed lines mark the ±2a envelope around the medians, where a is 
the root-mean-square scatter (~ 0.1 dex, due to software updates) estimated from the inter- 
quartile range. The median differences are larger than 0.1 dex only at the high-metallicity 
end i[Fe/H] > -0.6). 




Fig. A2. — Color-coded map in the top left panel shows the median difference between 
spectroscopic and revised photometric metallicities for ~ 50, 000 stars from SDSS Data 
Release 7. The median value is zero, and the root-mean-square scatter is 0.07 dex. The 
contours show the distribution of stars with r < 20 and at high Galactic latitudes. The top 
right panel shows the root-mean-square scatter of the difference between spectroscopic and 
photometric metallicities in each pixel. The top two panels are analogous to the bottom 
two panels in figure 2 from 108. The bottom left panel shows the photometric metallicity 
as a function of the spectroscopic metallicity. Individual stars are shown by small dots, and 
the median values of the difference are shown by large circles. The distribution of stars 
is shown as linearly-spaced contours. Note that the photometric metallicity saturates at 
[Fe/ H] ~ —2 at the low- metallicity end. The histogram in the bottom right panel shows the 
distribution of the difference between spectroscopic and photometric metallicities for stars 
with spectroscopic metallicity [Fe/H] > —2.2. A best-fit Gaussian centered on zero and 
with a width of 0.26 dex is shown by the dashed line. 
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Fig. A3. — Conditional metallicity distribution as a function of distance from the Galactic 
plane (top left), analogous to the bottom left panel of figure 9 from 108. Note that the new 
photometric metallicities include [Fe/H] > —0.4. The dashed line, which shows the median 
disk metallicity, is also revised (see text). The top right panel is analogous to figure 10 
from 108, and shows the metallicity distribution for stars with 5 < Z < 7 kpc, where Z 
is the distance from the Galactic plane. Note that the photometric-metallicity artifact at 
[Fe/H] = —0.5 discussed by 108 is no longer present. However, there is still evidence that 
disk stars exist at such large distances from the plane. The bottom two panels show the 
heliocentric rotational velocity for disk stars in two thin Z slices, and are analogous to the 
bottom right panel in figure 16 from 108. Note that the correlation between velocity and 
metallicity is still absent. 
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Fig. A4. — Analogous to figure 7 from 108. The symbols with error bars show the metallicity 
distribution for stars with 0.2 < g — r < 0.4, 7 < i? < 9 kpc, and distances from the 
Galactic plane as marked, where R is the Galactocentric cylindrical radius. The behavior 
is qualitatively similar to that seen in 108. The only significant quantitative difference is in 
the model for the metallicity distribution of disk stars (see text). 



